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SUMMARY 

Using  a  North  American  NAVION  as  a  flight  test  vehicle,  and  an 
electronic  analog  computer  as  a  simulator,  a  study  was  made  of  the 
feasibility  of  obtaining  certain  lateral  stability  derivatives  by  matching 
the  transient  responses  of  an  analog  computer  on  which  the  lateral 
operational  equations  of  motion  for  a  NAVION  aircraft  were  set  up,  to 
the  actual  airplane  lateral  response  in  sideslip  angle,  roll  rate,  and 
yaw  rate. 

The  matching  process  was  achieved  by  supplying  time  histories 
of  the  actual  recorded  control  motions  as  forcing  function  inputs  to  the 
computer,   and  systematically  varying  those  coefficients  of  the  analog 
which  corresponded  to  the  stability  derivatives  of  the  dynamic  equations 
of  motion. 

The  method  was   demonstrated  to  be  a  practical  procedure   for 
obtaining  the  lateral  stability  derivatives  of  a  linear  aircraft.      The   use 
of  four   different  types  of  control   deflections  as   forcing  functions   pro- 
vided a  means  for  matching  of  computer  responses  which  corresponded 
to  a  wide   range   of  flight  maneuver  s„ 
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The  method  is  considered  suitable  for   further  investigation,   and 
development  for  use   under   different  flight   conditions   and  with  non- 
linear  aircraft. 


IV 


TABLE  OF   CONTENTS 

Page 

Summary  iii 

Table  of  Symbols  ix 

Introduction  1 

Equipment  and  Procedure  5 

Equipment  5 

Airborne  Instrumentation  5 

Ground  Station  Instrumentation  1  3 

Instrument  Calibration  14 

Procedure  18 

Airborne  Phase  1 8 

Ground  Station  22 

Data  Reduction  23 

Discussion  26 

General  Description  of  Overall  Program  26 

Preliminary  Considerations  27 

Theoretical  Background  27 

Practical  Aspects  30 

Analog  of  the  Test  Aircraft  32 

Theoretical  Approximation  of  Stability  Derivatives  32 

Lateral  Equations  of  Motion  32 

Development  of  Method  of  Analysis  34 

Preliminary  Test  of  Feasibility  of  Method  36 

Application  of  Method  to  Flight  Test  Data  37 


TABLE  OF   CONTENTS   (Continued) 

Page 

Results  46 

Conclusions  and  Recommendations  50 

References  52 

Tables 

Figures 

Appendix  A     Dynamic  Response  of  Sideslip  Vane  A-l 

Appendix  B      Determination  of  Lateral  Stability  Deriva- 
tives  for   the   NAVION  from  Theoretical 
Formulae  and  from  Wind  Tunnel  and  Flight 
Tests   of  Geometrically  Similar   Airplanes  B-l 

Appendix  C      Corrections  to  Recorded  Aileron  Data  C-l 


vi 


TABLES 


I.  Methods  for  Analysis  of  Dynamic  Flight  Test  Data. 

II.  Physical  Characteristics  of  Test  Aircraft. 

Ill-a.         Effect  of  Variation  of  Stability  Derivatives  on 
Aircraft  Response. 

HI-D*         Stability  Derivatives  Which  Influence  Specific 
Aspects  of  Aircraft  Response. 

IV.  Primary  Effects  of  Individual  Stability  Derivatives 

on  Aircraft  Response. 

V.  Comparison  of  Theoretical  and  Final  Experimental 

Lateral  Stability  Derivatives  for  the   NAVION 
Aircraft. 


vii 


FIGURES 

1.  Three -View  Drawing  of  Test  Aircraft. 

2.  Photograph  of  Test  Aircraft. 

3.  Circuit  Diagram  for  Dynamic  Flight  Test  Instrumentation 
System  of  Test  Aircraft. 

4.  Photograph  of  Sideslip  Vane. 

5.  Instrumentation  System  Signal  Circuit  Schematic. 

6.  Photograph  of  Aircraft  Instrumentation. 

7.  Elements  of  a  Pulse  Width  Data  Collection  System. 

8.  Photograph  of  Ground  Station  Equipment. 

9.  Calibration  of  Roll  and  Yaw  Rate  Gyros. 

10.  Static  Calibration  of  Sideslip  Vane  Potentiometer. 

11.  Analog  of  NAVION  Lateral  Equations  of  Motion. 

12.  Potentiometer  Settings  for  Analog  of  NAVION. 

13.  Computer  Forcing  Function  Circuits. 

14.  Initial  Response  to  Rudder  Pulse,   Theoretical  Derivatives. 

15.  Matched  Response  to  Rudder  Pulse,    Final  Configuration. 

16.  Matched  Response  to  Rudder  Doublet,   Final  Configuration. 

17.  Matched  Response  to  Aileron  Doublet,    Final  Configuration. 

18.  Matched  Response  to  Aileron  Step  Function,   Final  Configuration, 
A-l.  Dynamic  Response  of  Sideslip  Vane. 


viii 


SYMBOLS 
CAPITAL  LETTERS 


b2/ 

AR  Aspect  Ratio,         'S 

CG  Center  of  Gravity 

CL  Lift  Coefficient,       Lift/^B 

Op  Drag  Coefficient,  '^s 

Cm  Pitching  Moment  Coefficient,    Pitching  Moment/   8 

Cg  Rolling  Moment   Coefficient,     Rolling  Moment/     , 

C  Yawing   Moment    Coefficient,     Yawing  Moment/^^ 

Cv  Lateral   Force   Coefficient,     Lateral  Force/a  a 

y  '^  ■ 


Wi«. 

_                   _____ 

cm 

_>  cy. 

% 

2>  (=> 

% 

% 

-^  >p 

c)    ___ 

cn 

'**«. 


Co,        _.      ^ 

3C, 


c 


*  **        S  %*, 


'-<vt 


C_.      _  __i 


3W 


c 

A  s,  -5 


IX 


"5C 


c^ 

V 

er          c3   '—  T\ 

D  Differential  Operator,       '  d     '*£- 

Also,  propeller  diameter,   ft. 

Ix  Moment  of  Inertia  about     X     axis,   slug  ft. 

I  a  Moment  of  Inertia  about     Z     axis,   slug  ft. 


Jx  Non- Dimensional  Inertia  Factor,        n  *x 

mb 


^— 2 


JH  Non -Dimensional  Inertia  Factor,      °)     *» 

b2m 

P  Period  of  Lateral  Oscillation,   seconds 

S  Wing  Area,  ft.2 

Sx  Characteristic  Area    (See  appropriate  subscript) 

Ti  Time  to  damp  to  one -half  amplitude,   sec. 


Tc  Thrust  Coefficient 


Vj[  Indicated  airspeed,  miles  per  hour 

V  True  airspeed,     "/sec. 

W  Aircraft  Test  Weight 


LOWER   CASE  LETTERS 

a  Lift   Curve   Slope,  per   radian 

b  Wing  Span,  ft 

C  Wing    Chord,   ft 

O  Mean  Aerodynamic  Chord,   ft 

e  Airplane    Efficiency   Factor 

ft/        2 
g  Acceleration  Due  to  Gravity,        '  sec 

h<j  Density  Altitude,  ft 

h  Pressure   Altitude,   ft 

Q  ^  Vertical  Tail  Length;    distance  from  center  of  pressure 

of  vertical  tail  to   CG,   ft 

X  sv  Sideslip   Vane   Length;    distance  from  sideslip  nave 

pivot  to  aircraft   CG,  ft 

m  Mass   of  Airplane,    slugs 

p  Rolling  Velocity,   radians/sec 

^C  Dynamic  Pressure,         '  Z^   ,  ft 

r  Yawing   Velocity,     rad/gec 

t  Time,   seconds 

Z^  Vertical   Tail  Height;    distance  from  center   of  pressure 

of  Vertical  Tail  to    X    axis 


XI 


<        •  • 


SUBSCRIPTS 


a 

Aileron 

fus 

Fuselage 

P 

Propeller 

r 

Rudder 

vt 

Vertical  Tail 

w  Wing 


%    x  Denote  First  and  Second  Derivatives  with  Respect 

i  to  Real  Time 


xii 


GREEK  LETTERS 

Oc.  Angle  of  Attack  of  Longitudinal  Body  Axis,   deg. 

r*  Angle  of  Sideslip,   radians 

"tf  Flight  Path  Angle,   deg. 

Control  Surface   Deflection 

Damping  Ratio 


S 


t^  Vertical  Tail   Efficiency  Factor 

\ 

A  Taper  Ratio,   Tip   Chord/Root  Chord 

^  Relative  Density  Factor,     m/psb 

\  Mass   Density  of  Air 

^  Time   Conversion  Factor,     m'osV 

S^  Angle   of  Bank,    radians 

4^  Angle   of  Yaw,    radians 

CO  Angular   Frequency,    radians/, 

^n  Natural   Frequency,     radians/, 

I  Dihedral  Angle,    degrees 

[\  Sweepback  of  Wing  Quarter-Chord  Line,   degrees 


sec 


sec 


xiii 


INTRODUCTION 

Since  the  close  of  World  War  II,   and  the  advent  of  the  high  speed 
and  altitude  characteristics  of  jet  aircraft,  the  need  for  a  quantitative 
knowledge  of  the  dynamics  of  flight  has  resulted  in  rapid  development 
of  the   field  of  aircraft  dynamics.     The  past   six  or   seven  years,    in 
particular,  have  brought  considerable  progress  in  the  field  of  dynamic 
stability  flight  testing  of  aircraft  and  missiles.     The  early  theoretical 
and  experimental  investigations  in  this  area  by  the   Cornell  Aeronau- 
tical Laboratory,   and  by  the  Instrumentation  Laboratory  of  the  Massa- 
chusetts  Institute   of  Technology,   have  been  supplemented  on    many 
fronts,   and  in  particular,  by  the   NACA,  which  in  recent  years  has 
developed  its   dynamic  stability  testing  techniques  to  the  point  where 
free-flight  testing  of  rocket  powered  and  telemeter  equipped  models 
has  become   commonplace. 

The  role  of  full  scale  dynamic  flight  testing  is  twofold.     First, 
it  provides  a  designer  with  the  conventional  data  on  the  characteristics 
of  an  aircraft  for  the  particular  conditions  of  the  test.    Secondly,  and 
in  the  long  view,  more  importantly,   dynamic  flight  testing  provides  a 
means  for  solution  of  the     inverse  problem     of  dynamics,   wherein 
the  equations  describing  a  system  are  assumed  known,  its  output  in 
response  to  suitable  forcing  functions  is  measured,   and  thus  the   co- 
efficients of  the  equations,   or  stability  derivatives,   can  be  determined. 


Once  established,   these  stability  derivatives  are  of  considerable  value. 
They  are  indispensable  to  the  aerodynamicist   who  must  calculate    sta- 
bility characteristics  over  a  wide  range  of  flight  conditions,   and  who  is 
thus   provided  with  a  measure   of  the   accuracy  of  the   theoretical  and 
empirical  derivatives  used  in  preliminary  design  calculations.     Deter- 
mination of  the    stability  derivatives   from  flight  tests   also  aids   the 
designer  by  indicating  lines  along  which  the  stability  may  be  improved, 
either  by  geometric  design  changes,   or  by  artificial  augmentation.     In 
the  latter  regard,  the  successful  design  of  any  system  for  automatic 
control  of  an  aircraft  is  greatly  dependent  upon  an  accurate  knowledge 
of  the    stability  derivatives    of  the   airplane. 

The   basic  approaches  to  the   problem  of  the   determination  of 
stability  derivatives  from  dynamic   flight  testing,   and  in  particular, 
from  experimental  time  histories,   can  be  separated  into  two  categories 
-  the  "equations  of  motion"  method,   which  essentially  solves  the   in- 
verse problem  described  above,   and  a  second  category,   "response 
curve  fitting,"  which  identifies  a  collection  of  methods  using  various 
mathematical  and  graphical  concepts  to  extract  from  the  data  certain 
combinations  of  stability  derivatives  and  inertia  terms  called  transfer 
function  coefficients .     All  of  the  individual  methods  which  make  up  the 
above  two  categories  will  not  be  described  in  this  report,   but  an  excel- 
lent discussion  of  each  current  method  can  be  found  in  Ref.  1.     Table  I 
presents  a  summary  and  comparison  of  these  various  methods    for 


analysis  of  dynamic  systems.     As  can  be  seen  from  Table  I,   the  choice 
of  a  method  for  analyzing  dynamic  response  data  will  depend  primarily 
on  whether  transfer  functions,   or  actual  stability  derivatives  are 
desired. 

Several  investigations  have  been  reported  (ie,  Ref.   2)  which 
have  attempted  to  obtain  stability  derivatives  by  applying  a  particular 
combination  of  the  above  referenced  methods  which  essentially  involved 
the  conversion  of  time  response  to  frequency  response,  followed  by  a 
matching  process  which  utilized  both  dynamic  response  templates  and 
an  analog  computer  to  produce  transfer  function  coefficients       These 
constants  were  then  shown  to  uniquely  determine  the  time  response   of 
an  analog  of  the  test  airplane  which  accurately  duplicated  the  time  his- 
tory of  the  actual  test  vehicle.     Evaluation  of  the  actual  aerodynamic 
stability  derivatives,  however,   required  supplemental  calculations,   and 
was  possible  only  in  the  case  of  those  derivatives  which  exerted  the 
strongest  influence  on  the  airplane  response. 

The  major  limitations  of  the  above  method  are  the  length  of  time 
required  for  conversion  of  the  flight  data  from  time  response  to  fre- 
quency response  form,   and  the  necessity  for  the  subsequent  use   of  tem- 
plates to  determine  approximate  transfer  function  coefficients.     It  was 
considered  that  both  of  these  limitations  could  be  avoided  in  the  applica- 
tion of  the  above  method  if  the  matching  process  were  to  be   applied 
directly  to  the  transient  flight  data. 


This  report,  then  is  concerned  with  a  study  of  the  feasibility  of 
obtaining  certain  lateral  stability  derivatives  by  matching  the  transient 
responses  from  an  analog  computer  on  which  the  operational  equations 
for  a  NAVION  aircraft  were  set  up,   to  the  actual  airplane  lateral   re- 
sponse in  sideslip  angle,   roll  rate,   and  yaw  rate.     The  matching  pro- 
cess was  achieved  by  supplying  time  histories  of  the  actual  recorded 
control  motions  as  inputs  to  the  computer,   and  systematically  varying 
those  coefficients  of  the  analog  which  corresponded  to  the  stability 
derivatives  of  the  dynamic  equations  of  motion. 

An  additional  objective  of  this  investigation  was  to  determine  the 
suitability  of  a  pulse  width  telemeter  system  as  a  data  collecting 
scheme  for  a  dynamic  flight  test  program  requiring  very  accurate 
determination  of  the  recorded  quantities. 

This  investigation  was  carried  out  at  the  James  Forrestal  Research 
Center,  Princeton  University,  Princeton,  New  Jersey,  during  the  1956-57 
academic  year. 


EQUIPMENT   AND  PROCEDURE 

EQUIPMENT 

Airborne    Instrumentation 

Test  Airplane.      The  test  airplane  was  a  North  American  NAVION, 
which  is  an  all  metal,   four  place,   low  wing,    single  engine  airplane. 
The  NAVION  is  equipped  with  a  Continental   E-185  engine,   rated  at  185 
HP  at  2300   RPM  and  29   inches  manifold  pressure,   at  sea  level.     The 
wing  is  unswept,   and  the  tail  configuration  conventionale     The   control 
surfaces  are  all  of  beaded  skin  construction,   and  have  no  trailing   edge 
extensions.     Aileron  and  rudder  trim  tabs  are  of  the  fixed  bend  type. 
The  elevator  tab  is  of  conventional  design,   and  can  be  adjusted  from  the 
cockpit. 

The  test  airplane  was  modified  by  the  removal  of  the  two   rear 
cockpit  seats  to  provide  space  for  the  test  instrumentation.     (Fig.    6). 
The  test  airplane  had  a  take-off  gross  weight  of  2870  pounds  with  full 
service,  pilot  and  co-pilot,   and  all  test  equipment. 

A  three  view  drawing  and  photograph  of  the  NAVION  are   shown 
as  Fig.  1  and  Fig.  2S   and  the  general  specifications  of  the  test  airplane 
are  tabulated  in  Table  II . 

Power  Supply.      Power  available  for  instrumentation  in  the   test 
airplane  included  both  12   and  28  volt   direct  current,    supplied  by  a  28 
volt  engine  driven  generator,   and  two  12   volt  batteries  connected  in 


series.     A   500  volt -ampere  inverter,   with  both  frequency  and  voltage 
regulation,   was  installed  behind  the  co-pilot's  seat  to  provide   115  volt 
400   cycle  three  phase  alternating  current  power  to  the  rate  gyros  and 
the  telemeter  transmitter. 

Considerable  effort  was  expended  in  obtaining  an  inverter  which 
could  be  frequency  controlled,    since  the  accuracy  of  the  rate  gyros  was 
dependent  upon  the  input  frequency  being  the  same  as  the  calibration 
frequency  of  400  cycles  per  second. 

A  circuit  diagram  of  the  aircraft   28  volt  system  and  the  inverter 
system  is  included  in  Fig„   3. 

Aircraft   Instrumentation.      Standard  aircraft  instruments  were 
used  to  determine  indicated  airspeed  and  altitude.     Measurement  of  the 
following  quantities  was  required  with  respect  to  stability  axes:     side- 
slip angle,   yaw  rate,   roll  rate,   and  aileron  and  rudder  position  angles. 
In  choosing  the  type  of  components  needed  to  measure  the  above  quan- 
tities,  an  attempt  was  made  to  utilize  instrumentation  compatible  with 
the  telemetering  equipment.     For  this  reason,   the  transducers  selected 
were  of  the  potentiometer  type,  not  only  because  they  were  capable  of 
providing  signals  of  a  sufficient  magnitude,  but  also  because  their  excita- 
tion source  readily  provided  zero  and  full  scale  reference   voltages. 

In  view  of  some  past  history  in  this  type  of  flight  testing  of  various 
amounts  of  coupling  between  the  response  of  the  aircraft  and  certain  of 
the  instrumentation  components,   considerable  attention  was  given  to  the 


choice  of  transducers  with  satisfactory  frequency  and  phase  response 
characteristics.  These  remarks  particularly  apply  to  the  choice  of  a 
suitable  device  for  recording  sideslip  angle. 

In  this  investigation,   sideslip  angle  was  sensed  by  a  light  vane 
mounted  on  a  boom  approximately  one  maximum  fuselage  diameter 
(four   feet)  ahead  of  the  left  wing  tip.     The  moment  of  inertia  of  the  vane 
was  very  low,   and  no  gearing  was  introduced  in  the  linkage  in  order  to 
keep  the  natural  frequency  and  damping  within  acceptable  limits.      An 
Electro-Mechanical  Laboratory  ultra-low  torque  potentiometer,   of  only 
.008    inch-ounces  friction  torque,   was  obtained  for  use  with  this   side- 
slip vane.     In  Appendix  A  is  presented  a  detailed  analysis  of  both  the 
theoretical  and  experimentally  obtained  natural  frequency  and  damping 
ratio  of  the  sideslip  vane. 

The  sideslip  vane  had  originally  been  designed  to  operate  through 
an  autosyn  transmitter,   and  thus  required  some  modification  in  order  to 
accommodate  a  potentiometer. 

hi  order  to  reduce  boom  resonance,   as  had  been  experience   in 
previous  tests  with  the  NAVION,  the  boom  for  the  sideslip  vane   was 
stiffened  by  three  angle  strips  equally  spaced  around  the  periphery,   and 
running  the  full  length  of  the  boom.     The  sideslip  vane  and  boom  are 
shown  in  the  assembled  condition  in  Fig.  4. 

Rudder  and  aileron  angles  were  measured  by  10,  000  ohm  Model 
5301  single  turn  Helipot  precision  potentiometers.     The  rudder  pot 


was  attached  to  the  rudder  itself,   and  was  actuated  by  a  cable  under 
spring  tension,,     At  the  time  the  instrumentation  system  was  designed, 
the  flight  test  program  was  planned  to  include  control  inputs  from  the 
rudder  onlyo     For  this  reason,   the  instrumentation  installed  for   re- 
cording aileron  angle  was  intended  only  to  indicate  if  the  ailerons 
remained  locked  throughout  the  rudder  input  and  the  resulting  tran- 
sient oscillation,,     Accordingly,  the  pickoff  of  aileron  angle  was  made 
from  control  cable  movement,   with  the  pot  attached  close  to  the  right 
aileron  in  order  to  avoid  the  effects  of  cable  stretch.     As  is  discussed 
in  Appendix   C,   the  later  decision  to  include  aileron  deflections   as 
forcing  functions  required  subsequent  correction  of  the  recorded  aileron 
data* 

The  measurement  of  yaw  and  roll  rates  was  made  by  type    36128 
Gianni  ni  rate  gyros.     This  type  of  gyro  operates  on  a  power  supply  of 
115  volts  400  cycle  three  phase  AC  to  the  motor,  and  115  volt   single 
phase   AC  to  the  heater.     The  transducer  is  of  the  potentiometer  type, 
with  a  DC  resistance  of  5250  ohms.     The  gyros  were  installed  as  near 
the  aircraft  center  of  gravity  as  possible.     However,   since  the   gyros 
were  of  the  angular  velocity  rather  than  acceleration  type,  there   was 
not  the  necessity  for  precise  location  that  is  associated  with  the  latter 
type  gyro.     A  level  and  inclinometer  were  used  to  align  the  gyros   so 
that  the  input  axes  would  be  parallel  to  the  appropriate  body  axis  of  the 
aircraft.     Since  deflection  of  the  gyro  spring -re  strained  gimbal  was 


only  two  degrees  at  full  range,  cross -coupling  effects  were  considered 
to  be  negligible.  However,  in  order  to  minimize  any  error  due  to  rate 
coupling,  the  spin  axis  of  each  gyro  was  aligned  parallel  to  the  body 
Y-axis.  The  gyros  were  bolted  directly  to  a  rigid  fuselage  member. 
Since  the  engine  vibration  frequency  was  sufficiently  above  the  natural 
frequency  of  the  gyros  (17  cycles  per  second),  there  was  no  engine 
noise  introduced  into  the  gyro  signal.  These  gyros  were  obtained  on 
loan  from  the  Instrumentation  Laboratory  of  the  U.  S.  Naval  Air  Test 
Center,   NAS  Patuxent  River,   Maryland. 

Signal    Circuit   Design*      The  component  of  the  instrumentation 
which  controlled  the  design  of  the  signal  circuit  was  the  sideslip  vane. 
The  requirements   for  recording  the  sideslip  angle   without  introducing 
phase  differences  were  negligible  torque,   and  a  relatively  high  natural 
frequency  of  the  vane.     The  transducer  selected  has  been  shown  to  have 
met  these  specifications   (Appendix  A),  but  the   fact  that  it  was   a   360 
degree  single  turn  potentiometer  pickup,   whereas  the  sideslip  vane  was 
expected  to  operate  through  only  a  moderate  arc,  necessitated  the  use 
of  45  volt   DC  excitation  of  the  sideslip  pot  rather  than  the  zero  to  five 
volts  range  supplied  to  the  remaining  transducers.    Since  the  resolution 
of  the  particular  sideslip  pot  used  was  exceptionally  good,  the  accuracy 
of  the    sideslip   system  was   comparable   to  that  of  a  system  utilizing 
autosyn  pickups,  but  had  none  of  the  problems  of  AC  signal   circuits 
and  was   completely  compatible   with  the  telemetering  and  recording 
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equipment.     The  use  of  a  45  volt   DC   signal  circuit  supply,   however, 
necessitated  a  method  of  supplying  this  voltage  to  the  sideslip  circuit, 
and  at  the  same  time,   limiting  the  full  scale  voltage  in  the  remaining 
signal  circuits  to  five  volts. 

In  achieving  the  above,   it  was  desirable  to  utilize  a  common  DC 
voltage  source  for  all  transducer  pickups,   since  the  telemeter  was  to 
sample  the  particular  supply  voltage  available   20  times  per  second, 
and  correct  this  to  full  scale  by  a  continuous  calibration. 

In  addition  to  the  above  requirements  on  the  signal  circuit  design, 
it  was  found  necessary  to  bias  the  sideslip  voltage  pot  in  order  that  the 
sideslip  signal  voltage  would  also  vary  between  zero  and  a  full  scale 
value  of  approximately  five  volts* 

The  particular  signal  circuit  design  employed  is  shown  in  Fig.   5, 

in  schematic  form.     The  essential  element  of  the  circuit  is  the   bias 

supply  voltage  potentiometer,  pot  "A".     With  the  sideslip  vane  against 

the  limit  stop  in  one  direction,  pot  "A"  was  adjusted  so  that  the  output 

voltage  of  the  sideslip  pot  was  zero.     The  vane  was  then  deflected  to  the 

opposite  limit  stop,  and  the  output  sideslip  pot  voltage  recorded.     Pot 

"B"  was  then  adjusted  to  bring  the  full  scale  reference  voltage,   and 

thus  also  the  voltage  in  the  remaining  signal  circuits,   in  agreement  with 

the  full  scale  sideslip  voltage.     This  procedure  was  repeated  two  or 

three  times,   as  necessary,   in  order  to  nullify  any  loading  on  the  system. 

Once  adjusted,   however,  pots    "A"  and   "B"   were  locked,   and  required 
no  further  adjustment  for  the  remainder  of  the  project. 
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During  the  initial  flight  and  transmission  tests,   considerable    AC 
ripple  was  found  to  enter  the   DC   signal  voltage  circuits,   very  likely 
from  the  400  cycle  inverter.     Condensers      C     and     D     in  Fig.    5 
were  installed,   and  were  found  to  reduce  the  ripple  to  a  negligible  value 
without  otherwise  altering  the  form  or  amplitude  of  the  signals . 

Fig.  6   shows  the  general  location  and  arrangement  of  the  air- 
borne instrumentation  components. 

Telemeter   System.      The  telemeter  unit  was  of  the  pulse  width 
type,   and  was  designed  and  built  by  the  Applied  Science  Corporation  of 
Princeton  (ASCOP).     This  unit  was  the  key  element  in  the  instrumenta- 
of  this  project.     The  basic  functional  components  of  this  pulse   width 
system  are  transducers,   all  of  the  potentiometer  type;   a  time  division 
multiplexer,   or  commutator;   and  electronic  coder,   or  keyer  to  put  the 
intelligence  into  pulse  form;   an  RF  transmitter;  and  ground  station 
equipment  for  receiving,   decoding,   and  recording  the  data  in  usable 
form. 

The  airborne  component  of  the  system  is  designated  the  ASCOP 
D  Series   PW  Multicoder,  and  the  ground  station  equipment  the  ASCOP 
M  Series   PW  Ground  Station.     A  schematic  diagram  of  the   functions 
of  the  pulse  width  data  system,   as  adapted  for  use  with  this  particular 
project,  is  presented  in  Fig.  7. 

A  unique  feature  and  very  great  advantage  of  pulse  width  data  sys- 
tems is  that  due  to  the  multiplexing  or  time  sharing  operation,  all  data 


12 


channels  traverse  the  same  path  during  transmission.     Since  zero  and 
full  scale  references  are  transmitted  with  the  data,  this  results   in  a 
calibration  signal,   effective  for  all  channels,   being  transmitted  20 
times  each  second,,     These  signals  supply  a  continuous  two-point  cali- 
bration,  and  thus  allow  for  the  removal  of  almost  all  of  the    system 
errors,   with  the  exception  of  errors  introduced  in  the  instrument  pick- 
ups and  the  mechanical  switch,   and  they  make  possible  the  use  of  this 
telemeter  in  tests  of  long  duration  and  varying  environmental  conditions. 

The  output  voltages  from  the  aircraft  instrumentation  pickups, 
and  the  calibration  voltages  described  above,   are  sampled  in  sequence 
by  a  rotary  switch.     At  the  normal  sampling  rate  of  20   samples  per 
second  on  each  of  the   45  channels,  the  data  pulse  length  for   any  one 
channel  is  short  in  comparison  to  the  interval  from  sample  to  sample. 
This  permits  the  use  of  a  storage  circuit  to  create  a  nearly  continuous 
voltage  output.     The  output  of  the  sampling  switch  is  a  series   of  data 
pulses  whose  voltage  amplitude  is  proportional  to  the  amplitude  of  the 
physical  quantity  measured  by  the  corresponding  pickup.     These  samples 
are  supplied  to  an  electronic  coder  circuit  (keyer)   which  converts  the 
amplitude  pulses  into  pulses  having  time  durations  proportional  to  the 
measured  quantities.     These  pulses  are  used  to  produce  modulation  of 
the  transmitted  RF  carrier.     The   RF   signal  is  detected  at  the   ground 
station,   with  the  output  of  the  receiver  being  a  train  of  pulses  identical 
in  nature  to  the  output  of  the  converter  in  the  airborne   equipment. 
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These  pulses  are  separated  by  a  switch  running  synchronously  with 
the  switch  at  the  transmitter,   with  the  pulses  representing  individual 
data  channels  being  placed  on  separate  outputs.     The  pulses  are  then 
demodulated  to  yield  continuous  voltages  representing  the  intelligence 
transmitted. 

Ground  Station  Instrumentation 


The  ground  station  equipment  included  the   Telemeter   Ground 
Station  (described  above),   a  Magnetic   Tape   Recorder,   a  GEDA 
Analog  Computer,   and  suitable  recorder s.     These  components   are 
illustrated  in  Fig.   8„ 

Magnetic  Tape  Recorder.      An  Ampex  Model   309 C  dual  track 
three  speed  tape  recorder  was  an  integral  component  of  the   ground 
station  equipment,   and  provided  a  convenient  means  for   recording, 
storage,   and  playback  of  the  standard  pulse  width  modulated  data. 
Provisions  were  also  included  for  placing  voice  transmissions  directly 
on  the  magnetic  tape  from  either  or  both  the  aircraft  and  the    ground 
station. 

GEDA  Analog   Computer.      The  analog  computer   used  in  this 
investigation  was  a  Goodyear   Aircraft   Corporation  Model   L3  (GEDA) 
linear  electronic  differential  analyzer.     Twenty-four  automatically 
stabilized  DC  computing  amplifiers  were  available  with  open-loop   DC 
gain  of  greater  than   5  x  10 '  ,  and  negligible  drift.      The   computer 
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incorporated  an  automatic   error   indicator,   and  had  a   guaranteed 
accuracy  of  one  percent.     Provisions   were  available   for   accurately 
setting  computer  board  potentiometers  by  the  use  of  a  special  null 
indicator . 

Data  Recorder.      Qualitative  analysis  of  the  flight  data  during 
tape  playbacks  was  made  on  a  Brush  two  channel  recorder.     Once  the 
particular  data  run  was  selected  for  further  analysis,   a  Sanborn  Model 
154-100B  four  channel  recorder  was  used.     This  instrument  had  a  very 
high  natural  frequency  (42   cps),   and  the  low  frequency  response   was 
flat  to  zero  cps„ 

Instrument   Calibration 

Five  components  of  the  instrumentation  system  required  calibra- 
tion.    These  were  the  two  rate  gyros,  the  two  control  deflection  poten- 
tiometers, and  the  sideslip  vane  potentiometer.     The  aircraft  altitude 
and  airspeed  indicating  systems  were  not  calibrated,  as  previous  cali- 
brations indicated  that  the  errors  in  these  systems  were  negligible  for 
the  conditions  of  this  investigation. 

The  two  rate  gyros  were  calibrated  on  a  pendulum  which  consisted 
of  a   3/4  inch  diameter  aluminum  shaft  and  a  12   inch  diameter  instru- 
ment pan.     The  pendulum,  the  overall  length  of  which  was   121.75  inches, 
was  pinned  to  the  overhead  truss  structure  of  the  Forrestal  hangar. 
A  scale  was  mounted  on  a  small  table  in  such  a  way  that  displacement 
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of  the  pan  from  its  rest  position  could  be  measured.     The  gyros  were 
mounted  one  at  a  time  on  a  wooden  base  plate  which  was  clamped  to 
the  pendulum  pan,   and  their  potentiometer  outputs  were  connected  to 
a  Brush  Pen  Recorder.     The  pendulum  was  displaced  from  its   rest 
position  and  allowed  to  swing  freely,,     Because  of  friction  in  the  system, 
swing  amplitudes  were  measured  on  two  successive  swings  in  the  same 
direction  and  averaged.    The  two  maximum  voltages  recorded  during 
the  swing  were  also  averaged  and  the  maximum  rate  attained  during 
the  swing  was  determined  from  the  following  relation: 


where 


P      -     natural  period  of  the  pendulum,   determined  from  the 
Brush  Recorder  trace,    seconds. 

J^     z      length  of  the  pendulum,   inches. 

d      -     average  maximum  displacement  of  the  pan  from  rest 
position,   inches. 


The  pendulum  was  allowed  to  swing  until  once  again  at  rest  with 
readings  taken  as  described  above  over  the  entire  range  of  gyro  rates. 
The  angular  rate  gradient  of  both  gyros  was  found  to  be  linear  and  equal 
to   0.610  degrees  per  second,  per  percent  full  scale  potentiometer 
voltage.     The  gyro  calibration  curves  are  presented  as   Fig.  9.     It  was 
found  that  the  gyros  were  very  sensitive  to  the  frequency  of  the   400 
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cycle  supply  voltage  and  as  a  result,    great  care  was  taken  during 
calibration,   and  also  during  the  actual  flight  test  program,   to  main- 
tain the  frequency  of  the  supply  voltage  at  exactly  400   cycles  per 
second. 

The  control  deflection  potentiometers  were  calibrated  by  meas- 
uring the  resistance  between  the  output  terminals  over  the  full  range 
of  control  deflections.     The  actual  control  deflections  were  measured 
with  a  vernier  protractor  capable  of  accuracy  to   0.1  degree.       The 
rudder  deflection  gradient  was  found  to  be  linear  and  equal  to   0.405 
degrees  per  percent  full  scale  potentiometer  voltage.     The  aileron  pot 
calibration  was  based  on  total  aileron  angle,   and  was  performed  over 
the  entire  range  of  aileron  travel  in  order  to  determine  the  extent   of 
aileron  differential  for  the  NAVION  aircraft.     The  calibration  was 
then  converted  to  average  aileron  angle,   and  was  found  to  be  linear  at 
a  constant  value  of  0.1933   degrees  per  percent  full  scale,   throughout 
the  deflection  range  used  in  the  flight  test,   which  included   -10  to  10 
degrees.     The  amount  of  aileron  differential  gearing  in  this  range  was 
found  to  be  extremely  small. 

The  sideslip  vane  potentiometer  was  calibrated  in  the  same 
manner  as  the  control  deflection  potentiometers,   and  the  calibration 
curve  is  presented  in   Fig.  10.     The  sideslip  gradient  was  found  to  be 
linear  and  equal  to   0.468  degrees  per  per  sent  full  scale  potentiometer 
voltage. 
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Since   the  telemeter  equipment   corrected   its  output   to  percent 
full  scale  voltage  and  since  full  scale  voltage  on  all   the  signal  circuit 
potentiometers  varied  by  the  same  increment  as  the  reference  voltage 
for  small  fluctuations  of  the   reference,   it  was   unnecessary  to   cali- 
brate the  potentiometers   through  the   telemeter    equipment. 

Careful  consideration  was  given  to  the  possible  errors   which 
might  be  introduced  in  the  data  by  the  dynamics  of  the  recording   in- 
struments.    The  natural  frequency  of  the  rate  gyros  was  17   cycles 
per  second  and  their  damping  ratio  was  between    .45  and   „7.      The 
natural  frequency  and  damping  ratio  for  the  sideslip  vane  were  experi- 
mentally determined  in  flight  to  be  15.44   cps  and  .081  respectively. 
(See   Appendix  A.)      Since  the  natural  frequency  of  the  aircraft   was 
approximately  0.4  cycles  per  second,  the  attenuation  of  the  signals 
due  to  the  dynamics  of  the  instruments  was  effectively  aero. 
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PROCEDURE 


General 


Aircraft  transient  responses  in  the  lateral  modes  of  motion  were 
obtained  by  introducing  rudder  and  aileron  control  motions  of  the  pulse, 
doublet,   and  step-function  types,,     Rudder  pulses,   rudder  doublets,   and 
aileron  doublets   were  found  to  be  the   most  effective   in  exciting  the 
*  Dutch  Roll"   motion,   which  in  turn  produced  the   most   satisfactory 
transient  time  histories  for  the  investigation  in  question.     Aileron  step- 
functions  were  also  used,  however,  primarily  to  provide  the  time 
history  of  roll-rate  response  to  lateral  control  motion. 

Ideally,    sufficient  transient  response  data  for  all  types  of  forcing 
functions  could  be  obtained  in  one  flight  of  approximately  an  hour's 
duration.     However,   due   primarily  to  the  need  of  the  investigators  to 
gain  some  proficiency  with  the  instrumentation  and  technique  of  applying 
satisfactory  forcing  functions,   several  flights  were  made  in  the  process 
of  collecting  data.     It  is  worthy  of  note   that  the  basic  scheme  of  instru- 
mentation and  data  collection  proved  singularly  reliable  during  the  flight 
test  -data  collection  phase.     After  initial  installation  the  instrumentation 
required  only  a  very  minimum  of  modification  and  adjustment,   and  in 
each  instance  the  trouble  was  minor  in  nature  and  easily  remedied. 
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Aircraft 

PrefLight  preparation  consisted  of  checking  the  full  scale  voltage 
applied  to  the  sideslip  potentiometer  to  insure  that  it  matched  that  of 
the  other  signal  potentiometers,   and  to  insure  that  the  full  scale  signal 
voltage  did  not  exceed  the  4-6  volt  range  which  the  telemetering 
equipment  needed  for  automatic  scale  calibration.     The  fuel  tank  was 
checked  and  if  necessary  topped-off  to  insure  that  the  aircraft   weight 
at  the  beginning  of  each  flight  was  constant.     When  deemed  advisable, 
a  complete  ground  check  of  all  communication  and  signal  transmission 
circuits  was  made  prior  to  take-off.     After  completing  the   ground 
checks  all  test  equipment  and  the  aircraft  generator  were   secured  to 
insure  maximum  engine  power  for  take-off. 

During  the  climb  to  test  altitude,  the  generator  and  inverter  were 
energized  and  power  was  applied  to  the  rate  gyro  heaters  and  to  the  air- 
borne telemeter  heater  and  filament  circuits.     Upon  reaching  test  alti- 
tude, voice  contact  was  established  with  the  telemeter  ground  station 
and  all  instrumentation  circuits  were  energized.     Because  of  an  appar- 
ent deficiency  in  the  inverter  frequency  control,  the  inverter  frequency 
was  monitored  continuously  and  maintained  at  400   cps  by  varying  the 
DC  voltage  supplied  by  the  aircraft  generator  to  the  inverter.     Voltage 
control  was  obtained  by  means  of  an  instrument -panel  mounted  potenti- 
ometer which  was  tied  into  the  aircraft  voltage  regulator  circuit   (Oper- 
ation of  the  synchronous  motors  in  the  rate  gyros  at  exactly  400  cps  was 
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necessary  to  insure  that  the  output    signals  followed  the  calibration 
curves  obtained  at   that   frequency). 

All  flight  tests  were  conducted  at  a  density  altitude  of  5000  feet 
and  at  a  true  airspeed  of  129.2   mph   (120  mph  IAS  at   5000     density 
altitude).     A  large  plot  of  free  air  temperature  versus  pressure   alti- 
tude was  carried  in  the  cockpit  to  facilitate  adjusting  pressure  altitude 
(altimeter  set  at  29.92  in.  Hg)   with  temperature  to  obtain  and  main- 
tain  5000  feet  density  altitude. 

The  general  technique  employed  in  the  actual  test  runs  fits   the 
following  description.     A  condition  of  steady-state  flight  was  established 
at  test  altitude  and  airspeed  and  held  for  at  least   30   seconds  prior  to 
each  control  input.     This  period  of  steady-state  flight  insured  that  the 
aircraft  was  as  nearly  trimmed  and  established  in  steady  state   condi- 
tions as  possible,   and  provided  a  sufficient  period  of  recorded  steady 
state  signals  to  allow  for  circuit  zero -signal  biasing  when  the  data  was 
played  back  into  the  analog  computer  or  pen  recorder.     At  the  end  of 
the  steady  state  period  a  radio  transmission  identifying  the  type  of  con- 
trol displacement  to  be  made  was  sent  and  recorded  on  the   magnetic 
tape.     Simultaneously,   all  controls  were  locked  mechanically  in  the 
cockpit  with  the  exception  of  the  control  to  be  displaced.     Five  seconds 
after  the  above  transmission  the  control  forcing  function  was  applied 
and  the  control  was  then  returned  to  the  original  steady- state  position 
and  locked.     All  controls  were  held  locked  until  the  resulting  aircraft 
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oscillation  had  essentially  subsided.     In  the  case  of  the  aileron  step 
functions  the  displacement  was  applied  as  rapidly  as  possible  and  the 
control  was  then  held  as   rigidly  as  possible.      The   resulting  motion 
in  this  case   was  permitted  to  proceed  for  4  to   6   seconds  before 
recovery  was  initiated. 

To  lock  the  elevator  and  ailerons  when  forcing  the  aircraft  with 
rudder  displacements,   a  chain  and  turnbuckle  system  attached  to  the 
deck  on  the  right  side  of  the  cockpit  was  used.     The  copilot's  control 
wheel  was  removed  and  replaced  by  a  sprocket  wheel  with  a  set  screw 
to  lock  it  to  the  control  shaft.     When  the  aircraft  was   established  in 
steady-state  flight,  the  sprocket  wheel  was  pushed  forward  on  the  con- 
trol shaft  until  it  butted  against  the  instrument  panel.     The  chain  was 
then  put  on  the  sprocket  and  the  turnbuckle s  taken  up,  making  the  chain 
taut.     Just  prior  to  initiating  the  rudder  deflection,  the  set  screw  was 
tightened,  locking  the  sprocket  fitting  to  the  control  shaft.     With  the 
chain  on  and  tight,  the  ailerons  were  effectively  locked,  and,   with  the 
sprocket  fitting  against  the  instrument  panel,  no  forward  motion  of 
the  control  shaft  was  possible.     The  copilot  then  held  forward  pressure 
on  the  sprocket  fitting  to  prevent  aft  movement  of  the  control  shaft. 

To  lock  the  rudder  when  putting  in  aileron  displacements,   and  to 
return  the  rudder  to  its  steady-state  position  after  rudder  displace- 
ments,  a  chain  and  turnbuckle  system  on  the  pilot's  right  rudder  pedal 
was  employed.     This  was  attached  to  the  right  hand  seat  track  on  the 
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pilot  fe  side.     Because  the  aircraft  required  a  slight  amount  of  left 
rudder  for  trim  at  the  test  flight  condition,  this  attachment  was  also 
useful  in  establishing  trimmed  flight  by  adjusting  the  turnbuckle. 

Ground  Station 

In  order  to  visually  monitor  the  data  runs,   an  external  patch- 
board provided  on  the  telemeter  receiver   (See  Fig.   8)  was  wired  in 
such  a  way  as  to  utilize  a  voltmeter -type  translator  for  each  of   the 
recorded  variables,  and  also  to  allow  the  decoded  output  to  be  placed 
on  a  pen  recorder.     The  telemeter  receiver  and  the  magnetic  tape 
recorder  had  been  integrated  prior  to  this  investigation,  thus  requiring 
only  the  proper  positioning  of  controls  on  the  telemeter  console  to 
accomplish  tape  recording  of  the  flight   data.     Each  quantity  measured 
during  the  flight  tests  was  transmitted  on  three  separate  telemeter 
channels.     Thus,  each  measured  variable  was  sampled  60  times  each 
second.     Once  the  airborne  equipment  was  operating,  the  telemeter 
receiver  was  tuned  for  maximum  signal  reception,  and  the     zero     and 
full-scale     reference  voltages  were  adjusted  for  proper  calibration  of 
output  signals. 

All  pertinent  communications  between  the  ground  station  and  the 
aircraft  were  recorded  on  magnetic  tape  along  with  the  flight  test  data. 
Prior  to  commencing  the  recording  of  test  data  on  each  flight,    the 
ground  station  operator  pre-recorded  on  the  magnetic  tape  the   data 
necessary  to  identify  the  particular  flight,   and  test  conditions  as 
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regarded  aircraft  weight,  lift  coefficient,   and  the  telemeter   channels 
being  used  for  each  measured  quantity. 

Concurrent  with  the  reception  and  tape   recording  of  the   flight 
test  data,   a  two  or  four  channel  pen  recorder  was  used  to  present 
visually  for  monitoring  all  or  a  part  of  the  test  data  as  the   maneuvers 
were  being  performed.     Thus,   if  the  control  input  did  not  produce  well 
excited  responses,   or  if  a  control  Burface  that  was  to  be  held  locked 
moved  inadvertently,   or  if  the  forcing  control  was  not  returned  to,   and 
held  at,  the  original  steady-state  position  this  was  immediately  appar- 
ent.    In  any  of  these  eventualities  the  run  could  be  rejected  immediately 
on  a  qualitative  basis  and  the  test  repeated. 

Data  Reduction 

The  only  flight  variable  requiring  correction  to  the  measured 
quantity  was  sideslip  angle   {f?  )   which  had  to  be  corrected  for  location 
of  the  sideslip  vane  ahead  of  the  aircraft  center  of  gravity.     The  required 
correction  was  : 

'     corr.  /     measured  W      '      ' 

In  order  to  have  all  flight  variable  signals  in  the  same  angle- 
time -voltage  relationship,  the  flight  data  was  processed  through  scaling 
circuits  on  the  analog  computer  prior  to  plotting  on  a  pen  recorder. 
The  scaling  factors  used  were  determined  from  the   calibration  curves 
for  each  of  the  signal  potentiometers.     Where  necessary,  the  signs   of 
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the  measured  variables  were  changed  to  correspond  to  the  sign  of  the 
corresponding  output  of  the  analog  probelm.     A  .25  uf  capacitor  was 
placed  across   each  of  the  output  terminals  of  the  telemeter.      This 
size  capacitor  was  found  to  help  in  eliminating  high  frequency  noise  in 
the  flight   data  signals  while  at  the  same  time  having  negligible  effect 
on  the   amplitude  and  phase   of  the   signals.     The   computer   circuits 
used  in  correcting  the  sideslip  angle  and  scaling  the  flight  test   data 
were  as  follows: 
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Flight  Data  Scaling  and  Correction  Circuits 


In  order  to  use  the  actual  flight  test  control  deflections  as  forcing 
functions  for  the  analog  problem  it  was  necessary  to  feed  them  into  the 
computer  through  a  bias  circuit  so  that  the  steady-state  control  position 
signals,   which  were  of  the  order  of  -50  volts  at  the  telemeter   output 
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terminals,   could  be  reduced  to  a  zero  voltage  input  to  the   analog 
problem.     A  potentiometer  was  also  required  in  each  input  circuit  to 
adjust  the  amplitude  of  the  control  deflection  voltage  to  match  the 
scaled  magnitude  of  the  control  deflection  on  the  flight  test   data  plots. 
As  in  the  data  reduction  circuits,   a  .25  uf  capacitor  was  put  across  the 
telemeter  output  terminals.     The  above  described  circuits  are  included 
in  Fig.  13. 
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DISCUSSION 

General  Description  of  Overall  Program 

In  obtaining  flight  data,  four  different  types  of  control  deflections 
were  used  to  excite  the  lateral  modes  of  the   NAVION.     These  were   a 
rudder  pulse,   a  rudder  doublet,   an  aileron  doublet,   and  an  aileron  step 
function.     By  means  of  telemeter  equipment,  the  transient  response  of 
the  aircraft  as  well  as  the  control  deflections  as  a  function  of  time  were 
recorded  on  magnetic  tape  at  the  ground  station.     From  a  large    selec- 
tion of  such  data,   one  representative  trace  was  made  of  each  of  the  four 
types  of  forcing  functions  mentioned  above.     This  was  accomplished  by 
processing  the  flight  data  through  scaling  circuits  in  the  analog  computer 
and  recording  the  variables  in  time  history  form  on  a  Sanborn  recorder. 
A  transparent   overlay  was   then  made   of  each  of  the   four   traces   of 
scaled  flight  data  in  order  to  facilitate  comparison  of  the  flight  and  com- 
puter  responses. 

Next,  the   lateral   stability  derivatives   of  the   NAVION  were 
theoretically  estimated.     The  lateral  equations  of  motion  based  on  these 
derivatives  were  placed  on  an  analog  computer  in  such  a  way  as  to  enable 
independent  variation  of  each  stability  derivative.     The   magnetically 
recorded  signal  corresponding  to  one  of  the  four  chosen  control  deflec- 
tions was  then  fed  directly  into  the  computer  to  force  the  analog.    Both 
the  character  and  magnitude  of  the  forcing  function  signal  were  checked 
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by  means  of  the  overlay  of  the  flight   data.     Once  the  forcing  function 
was  properly  adjusted,  the  computer  responses  to  the  forcing  function 
were  recorded  and  compared  to  the  flight  responses  by  means  of  the 
same  overlay.     Adjustments  were  then  made  to  the  computer   deriva- 
tives as  necessary  to  improve  the  match  between  the  flight  data  and 
the  computer  response „     The  procedure  of  forcing  the  computer   was 
then  repeated,   using  another  one  of  the  four  forcing  signals,   and  mak- 
ing appropriate  adjustments  to  derivative  values.     This  procedure 
proved  to  be  a  converging  iteration  process  when  repeated  several 
times,   using  all  four  types  of  forcing  functions.,     When  a  set  of  deriva- 
tives was  found  such  that  the  response  of  the  computer  gave  satisfactory 
matching  with  the  flight  data  for  all  four  of  the  forcing  functions,   and  no 
further  improvement  seemed  justifiable  in  relation  to  the  amount  of  time 
involved,  the  value  of  each  stability  derivative  was  read  out  of  the  com- 
puter and  tabulated. 

Preliminary    Considerations 

Theoretical  Background,,     At  the  onset,  of  the  investigation,    it 
was  necessary  to  make  a  purely  theoretical  analysis  of  the  proposed 
method  of  attack  of  the  problem.     In  effect,  the  pertinent  questions  were: 
a)    If  a  set  of  coefficients  is  determined  which,  when  placed  in  the 
equations  of  motion,  duplicate  the  response  of  the  aircraft, 
can  it  be  proved  that  these  coefficients  are  unique  rather  than 
just  one  set  of  an  entire  family  of  possible  solutions? 
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b)    With  two  possible  methods  of  exciting  the  lateral  mode,   i.e. , 
the  rudder  and  the  ailerons,  how  much  data  should  be  col- 
lected and  what  should  be  the  character  of  the  control  deflec- 
tions used  to  excite  the  motion? 

Examining  these  questions,   it  was  determined  that  one  was  related 
to  the  other.     The  answer  to  question  (a)   is  a  function  of  the   amount  of 
independent  data  available.     Should  there  be  a  case  where,   in  fact,   an 
entire  family  of  solutions  exists,   then  there  would  be  an  interdependence 
between  at  least  two  of  the  coefficients,   and  one  coefficient  could  be   ex- 
pressed in  terms  of  one  or  more  of  the  other  coefficients.     If  it  were 
then  possible  to  obtain  more  independent  data,   a  point  would  be  reached 
where  there  were  sufficient  conditions  to  limit  the  family  of  solutions  to 
a  single  solution.     Applying  the  latter  case  to  the  investigation  at  hand, 
it  would  then  be  known  that  if  a  set  of  coefficients  were  obtained  by  the 
response -matching  technique,  these  coefficients  would  be  unique,   i.e., 
the  only  set  of  coefficients  which  would  produce  the  matching  of  re- 
sponses.    The  stability  derivatives  contained  in  the  lumped  coefficients 
would  then  indeed  be  the  stability  derivatives  of  the  aircraft.      The 
amount  of  independent  data  available  would  be  in  part  a  function  of  the 
control  deflections  used,   and  in  addition  would  depend  on  the  number  of 
quantities  such  as  roll  rate,   yaw  rate,   and  sideslip  angle  which  were 
recorded  during  the  response. 

Fortunately,   a  theoretical  investigation  into  this  subject  was 
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available  in  Ref.  4.     The  authors,   Messrs.   Klein  and  Sedney„     assume 
that  the  equations  of  motion  are  linear  with  constant  coefficients    (how- 
ever,  it  is  pointed  out  that  the  same  method  of  analysis  can  be  applied 
to  more  general  types  of  equations),   and  also  that  the  data  is  free  from 
experimental  error.     The  results  of  their  analysis  is  a  simple  working 
rule  which  can  be  stated  in  words  as  follows:     In  each  equation  describ- 
ing a  system,  the  number  of  unknown  coefficients  which  can  be  deter- 
mined plus  the  number  of  known  coefficients  which  may  be  checked  are 
equal  to  the  total  number  of  measured  quantities  less  the   number   of 
equations . 

This  result  was  applied  to  the  lateral  equations  of  motion  of  the 
NAVION,   with   Jx,    Jz,    Jxa,   and  CL  known  .      The  quantities  which 
were  to  be  recorded  during  the  testing  were  roll  rate,   yaw  rate,   side- 
slip angle,   and  control  deflections       In  application  of  the   Klein  and 
Sedney  criterion,   roll  angle,   roll  acceleration,   yaw  angle,   yaw  accel- 
eration,  and  sideslip  rate  are  considered  to  be  recorded  in  a  concealed 
form  in  the  other  time  histories,    since  they  are  obtainable  by  differen- 
tiation or  integration  of  the  appropriate  trace.     By  following  the  working 
rule  of  the  criterion,   it  was  found  that  in  each  of  the  three  equations  of 
motion  there  was  actually  an  excess  of  independent  data  under  the  re- 
cording conditions  proposed,   and  therefore  the  solution  for  the  coeffi- 
cients would  be  unique. 
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Practical   Aspects.      Although  no  theoretical  limitation  exists  in 

the  case  of  the  lateral  mode,  there  is  a  practical  limitation  which 

applies  to  the  analysis  of  all  dynamic  response  data,,     The   following 

statements  are  extracted  from   Chapter   10   of  Ref.  1: 

Instrument  errors  and  unknown  atmospheric  disturbances 
are  usually  present  to  a  sufficient  degree  to  introduce  significant 
errors  into  the  solution  of  the  inverse  problem.     As  a  general 
qualitative  rule,   the  accuracy  with  which  one  can  obtain  each  of 
the  stability  derivatives  is  associated  with  the  importance  of  that 
stability  derivative  in  influencing  the  motion  of  the  aircraft  dur- 
ing the  test,  the  highest  accuracy  being  obtained  for  the  more 
important  stability  derivatives. 

Consideration  of  the  means  for  improving  the  accuracy 
with  which  the  less  important  derivatives  can  be  obtained  leads 
to  these  three  approaches  to  the  problem; 

1.  Overall  improvement  of  the  accuracy  of  the  recording 
instrumentation  and  improvement  of  the  calibration 
and  pre -flight  testing  techniques. 

2.  Selection  of  test  maneuvers  which  increase  the  force 
and  moment  contributions  due  to  the  smaller  deriva- 
tives compared  to  the  forces  and  moments  due  to  the 
other  derivatives. 

3.  Performance  of  flight  tests  under  so-called  ideal 
conditions,  i.e.,   air  is  smooth  and  stable. 

Selection  of  flight  test  instrumentation  components  was  made 
with  the  requirements  of  (1)   above  in  mind.     As  previously  discussed 
under     'Airborne  Instrumentation,      the  transducers  were  chosen  to  be 
compatible  with  the  recording  instrumentation  in  order  to   eliminate 
excessive  electronic  modification  of  the  signal  voltages  prior  to  record- 
ing.    This  effectively  reduced  possible  sources  of  error.     Calibration 
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was  carefully  performed  and  the  resolution  and  linearity  of  the  com- 
ponents are  readily  apparent  from  the  calibration  curves,   as  presented 
in  Fig.  9   and  Fig.  10.     Test  flights  were  performed  under  conditions 
of  smooth  and  stable  air. 

The  selection  of  suitable  forcing  functions  for  the  lateral  mode 
was  initially  considered  from  the  viewpoint  of  Klein  and  Sedney,   who 
state  that  theoretically,   the  response  of  the  aircraft  to  a  single  forcing 
function  is  sufficient  to  establish  the  stability  derivatives,   with  the 
exception  of  the  control  derivatives .     Accordingly,   a  rudder  pulse  was 
selected  as  the  forcing  function  most  suitable  for  excitation  of  the 
lateral  mode.     However,   the  prospect  of  extending  the  method  of  analy- 
sis to  obtain  all  of  the  control   derivatives,   plus  the  possibility  of 
improving  the  accuracy  of  certain  derivatives  through  the  use  of  addi- 
tional forcing  functions,   was  made  very  attractive  when  considered 
from  the  standpoint  of  paragraph   (2)   of  the  above  quotation  from  Ref.  1. 
Four  forcing  functions  were  therefore  decided  upon.     Considering  the 
rudder  pulse  as  the  basic  data,   a  rudder  doublet  was  utilized  to  in- 
crease the  amplitudes  of  the  lateral  responses,   and  to  reduce   the 
tendency  of  the  aircraft  to  fall  off  on  one  wing  during  the  transient  as 
had  been  observed  in  the  case  of  excitation  by  a  rudder  pulse.      The 
aileron  doublet  was  considered  as  an  effective  means  of  evaluating 
CO      ,     Ch       ,    and  the  aileron  derivatives,   whereas  a  time-history  of 
the  roll  rate  response  to  an  aileron  step  function  was  considered  to  be 
the  best  means  for  analyzing  the  stability  of  the  spiral  mode. 
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Analog   of  the    Test  Aircraft 

Theoretical  Approximation  of  Derivatives.      In  order  to  construct 
an  analog  of  the  lateral  equations  of  motion,   it  was  necessary  to  esti- 
mate the  value  of  the  stability  derivatives  from  whence  the  approximate 
value  of  each  lumped  coefficient  could  be  determined.     Anticipating  the 
possibility  of  extremely  slow  convergence  of  the  iteration  process  used 
on  the  computer,  the  theoretical  value  for  each  derivative  was   deter- 
mined as  accurately  as  possible  by  several  different  methods  and  then 
a  median  value  was  taken.     The  results  and  a  complete  explanation  of 
this  work  is  contained  in  Appendix  B. 

Lateral  Equations  of  Motion.      The  lateral  equations  of  motion  of 
a  linear  aircraft  are  as  follows: 


=1 
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For  ease  of  construction  of  the  analogy  the  equations  can  be 
rewritten  incorporating  the  following  changes: 

a)  in  terms  of  the  highest   derivative. 

b)  in  real  time  by  the  explicit  insertion  of  the  non-dimensional 
time  constant ,  IT  . 
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c)     for  simplicity,   the  algebraic  signs  changed  where  necessary 
so  that  all  work  from  this  point  on  could  be  conducted  in 
terms  of  the  absolute  value  of  each  derivative., 

The  resulting  equations  are: 


x 
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The  numerical  values  of  the  non-dimensional  constants  and  the 
estimated  derivatives  were  then  substituted  into  the  equations „      The 
terms  containing  Jxz   were  found  to  be  so  small  in  comparison  to  the 
other  terms  that  they  were  immediately  discarded  as  negligible: 


<*- 


The  constructed  analog  for  the  above  equations  is  shown 
schematically  in  Fig.  11 « 
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Development   of  Method  of  Analysis 

Prior  to  the  analysis  of  the  actual  flight  test   data,   a  step  by 
step  procedure  was  developed  for  matching  the  computer  response  to 
flight   data.     First,   an  investigation  was  made  to  determine  the  effect 
of  each  stability  derivative  on  the  response  of  the  system  in  terms  of 
amplitude,  period,   and  damping.,     The  analog  of  the  test   airplane 
(Fig.  11)   was  placed  on  the  analog  computer,   and  the  potentiometers 
set  to  correspond  to  the   NAVION  theoretical  stability  derivatives   of 
Table   V.     An  electronic  function  generator  was  used  to  force  the 
analog  with  a  square  pulse  type  function,  to  approximate  the  form  of  a 
rudder  pulse.     A  series  of  runs  was  then  made  during  which  the  sta- 
bility derivatives  were  varied  successively  by  a  factor  of  two  and  then 
returned  to  their  original  theoretical  values.     From  these  tests,    the 
change  in  the  response  due  to  the  change  of  any  one  derivative   was  re- 
corded,  and  the  results  included  in  Table   ILt-a.     A  cross -table   was 
then  prepared  showing  which  derivatives  should  be  varied  to  obtain  a 
desired  change  in  the  response  of  the  system.     These  results   are 
shown  in  Table   Ill-b. 

With  the  aid  of  the  information  presented  in  Tables   ni-a  and 
m-b,   an  iterative  procedure  for  the  matching  of  computer  response  with 
flight   data  was  outlined  as  follows: 

(a)     Compare  flight   data  and  computer  response  curves  obtained 
from  the  same  forcing  function.     Note  discrepancies   in 
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amplitude,   phase,    damping,   etc.,   between  corresponding 
responses. 

(b)  Concentrate   on  the   major   discrepancy.     By  reference  to 
Tables   Ill-a  and  Ill-b,   determine   which  derivative,    or 
derivatives,   could  be   altered  to  correct  this  type   of 
discrepancy. 

(c)  Iteratively  apply  the   corrections   indicated  in  (b)   above, 
commencing  with  the  derivative  having  the  most  pronounced 
effect  on  the  particular   discrepancy.     When  no  further 
improvement  can  be  observed  by  adjusting  the  value  of  the 
first   derivative,   continue  in  a  consecutive  adjustment  of  all 
the  derivatives  affecting  the  discrepancy. 

(d)  Repeat   (a),    (b),   and  (c)   above  for  additional  discrepancies 
until  general  agreement  is  obtained  between  flight   data  and 
the  computer  response  curves. 

(e)  Continue  by  comparing  flight  data  and  computer   response 
curves  obtained  from  a  second  forcing  function.     Again, 
repeat   (a),    (b),   and  (c)  until  agreement  is  obtained. 

(f)  Repeat   (e)   for  as  many  forcing  functions  as  are  available  or 
practical.     When  the  response  to  all  forcing  functions  have 
been  examined  once,  return  to  the  first  forcing  function  and 
repeat  the  entire  process  outlined  above;   continue  until  a 
set  of  derivatives  is  obtained  which  give  the  best  possible 
response  match  to  all  forcing  functions  used. 
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Preliminary  Test   of  Feasibility  of  Method 

A  program  was  next  performed  to  establish  the  feasibility  of 
using  an  analog  computer  to  curve -fit  flight  data  and  yield  unique 
stability  derivatives.     It  was  considered  that  if  the  transient  response 
of  the  computer  to  a  rudder-pulse  type  of  forcing  function  were  record- 
ed,  and  considered  as  error -free     "flight   data"      and  the   computer 
potentiometers  corresponding  to  stability  derivatives  then  changed, 
that  subsequent  re -matching  of  the  original  synthetic   "flight   data" 
response  should  yield  computer  potentiometer  settings  identical  to 
those  of  the      flight  data.      It  was  found  that  this  type   of  analysis 
produced  excellent  results,   with  the  exception  of  the  determination  of 
^Sy-     '    **Ua  '    an(i   *~"y&       •      ^he  responses  were  found  to  be  too 

insensitive  to  changes  in  these  three  derivatives  in  order  to  establish 
them  with  satisfactory  accuracy.     The  final  value  of    C.,     contained 
an  error  of  approximately  ten  percent.     All  other  derivatives   were 
obtained  within  an  error  of  five  percent.       C^  and    Cy  were 

obtained  with  virtually  no  error  ;    this  latter  result  was  anticipated, 
from  consideration  of  the  large  relative  magnitude  of  the  lumped 
coefficients  containing  these  terms. 

Unfortunately,  the  abore  program  was  not  extended  to  the  use  of 
an  aileron  input.     Although  the  results  were  gratifying  for  the  rudder 
forcing  function  used,  they  might  well  have  been  improved  if  the  addi- 
tional forcing  function  data  had  been  utilized.     In  particular,  the  use  of 
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aileron  forcing  functions  might  have  resulted  in  a  better  indication  of 
the  effect  of  the  derivative    C  ^ 

The  investigators  considered  that  the  results  of  this  preliminary 
program  established  the  feasibility  of  the  basic  method.     It  was   also 
encouraging  to  note  that  with  a  minimum  of  practice,  the   iteration  pro- 
cedure developed  for  matching  the  analog  responses  to  flight   data 
resulted  in  moderately  rapid  convergence.     With  ten  stability  deriva- 
tives all  acting  as  variables  in  the  inverse  problem,   a  great   deal  of 
difficulty  had  been  anticipated  in  effecting  a  convergence  when  using 
what  is  basically  a  simple  iteration  process,   but  to  the  degree  of  accur- 
acy obtained  in  the  results  mentioned  above,   no  significant  problems 
were  realized,, 

In  future  applications  of  this  method,  the  procedures  described  in 

it  .  "  " 

the  sections      Development  of  Method  of  Analysis      and     Preliminary 

Test  of  Feasibility  of  Method     would  be  useful  to  an  investigator  from 

two  standpoints:    they  would  provide  data  on  the  sensitivity  of  the  specific 

configuration  under  investigation  to  changes  in  individual  derivatives, 

and  would  provide  valuable  experience  in  manipulating  the  coefficients 

to  obtain  the  desired  response. 

Application  of  Method  to  Flight   Test   Data 

The  four  control  deflections  which  were  selected  and  the  corres- 
ponding responses  of  the  aircraft  in  the  three  variables  /$  ,     W     ,  and 
y>      are  shown  in  Figs.  15-18.     As  mentioned  previously,   a  transparent 
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overlay  was  made  of  each  of  these  four  sets  of  data  to  facilitate   com- 
parison of  the  computer  response  with  the  flight   data. 

One  of  the  important  features  of  the  instrumentation  was  the  fact 
that  signals  proportional  to  the  four  control  deflections  were  available 
for  forcing  the  computer  by  playing  the  recording  of  flight  data,   stored 
on  the  magnetic  tape,  through  the  decoding  component  of  the  telemeter 
ground  station.     From  the  output  of  the  telemeter,  the  circuity  shown 
in  Fig.  13   was  used  to  bias  the  steady  state  signals  of  the  controls  to 
zero  and  to  provide  proper  scaling  of  the  control  deflection  signal 
magnitudes.     Without  this  feature,  it  would  have  been  necessary  to 
obtain  the  forcing  function  for  the  computer  from  an  electronic  function 
generator,   with  limited  prospects  for  exact  duplication  of  the  aircraft 
control  deflections.     The  use  of  the  exact  forcing  function  from  the  flight 
data  thus  eliminated  a  possible  source  of  error  by  maximizing  the  com- 
patibility of  the  data  collection  system  and  the  analog  computer. 

The  relationship  between  the  actual  aircraft  response  and  the 
response  of  the  analog  containing  the  initial  theoretical  derivatives, 
when  each  was  forced  by  a  rudder  pulse,  is  shown  in  Fig.  14.     With  this 
initial  analog  response   (Fig.  14)   as  a  starting  point,  the  attempt  to 
curve -fit  the  analog  to  the  flight  data  was  commenced. 

The  initial  steps  of  the  curve  fitting  process  were  the  same  as  the 
first  four  parts  of  the  procedure  previously  described  in     Development 
of  Method  of  Analysis.        Reference  to   Fig.  14  indicated  that  the  most 
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obvious  discrepancy  was  in  the  period  of  the   Dutch  roll  oscillation. 

Since   Tables   III-a  and  m-b   showed    C  ^  to  be  the  derivative  most 

effective  in  adjusting  the  period,  the  value   of    C  y>  A       was  increased 

by  small  increments  until  a  relatively  good  period  match  was  obtained. 

The  values  of  the  rudder  control  derivatives    (Cj.         ,     Co/       ,   C,?,       ) 

were  then  corrected  to  reflect  the  change  in    C  ^  .      The  next  steps 

P 

were  to  attack  successively  the  discrepancies  in  damping,    yaw-rate 
amplitude,   and  roll-rate  amplitude  by  adjusting  the  appropriate  deriva- 
tives in  accordance  with  the  information  in  Tables   ni-a  and  III-b. 
This  procedure  improved  the  match  between  the  analog  and  the  aircraft 
responses  appreciably.     However,  in  an  attempt  to  further  improve  the 
matching,   and  at  the  same  time  gain  further  insight  into  the  interplay 
of   the  effects  of  the  various  derivatives,   all  derivatives  not  adjusted 
in  the  foregoing  procedures  were  successively  altered  to  observe  their 
effect  on  the  analog  responses.     When  no  change,  or  an  adverse  change, 
in  the  responses  was  encountered,  the  derivative  under   investigation 
was  returned  to  its  original  value.     Those  changes  which  appeared 
beneficial  were  retained.     In  each  case,  the  effect  of  altering  a  deriva- 
tive,  and  the  extent  of  the  change  were  recorded. 

After  a  reasonably  close  match  between  the  analog  and  aircraft 
responses  to  a  rudder  pulse  was  obtained,  the  analog  was  next  forced 
by  a  rudder  doublet,   with  the  computer  potentiometers  still  set  on  the 
derivatives  established  by  matching  responses  to  the  rudder  pulse. 
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The  procedure  which  was  then  followed  to  minimize  the  discrepancies 
between  the  analog  and  aircraft  responses  was  identical  to  that  used 
for  response  to  a  rudder  pulse,   and  again  produced  an  improvement  in 
the  analog  responses. 

The  responses  to  an  aileron  doublet  and  aileron  step  function 
were  then  investigated  in  a  like  manner;  that  is,  by  first  concentrating 
on  the  derivative  changes  indicated  in  Tables  ni-a  and  Ill-b,  and  then 
investigating  the  individual  effects  of  all  other  derivatives.  During  the 
initial  analysis  of  the  aileron  data,  it  was  determined  that  an  error 
existed  in  the  aileron  signal  due  to  cable  stretch  in  the  control  system. 
Corrections  required  for  this  error  are  discussed  in  Appendix   C. 

During  the  conduct  of  the  actual  flight  tests,   it  had  been  noted 
that  in  response  to  an  aileron  step  function,  the  aircraft   developed  a 
nose-down  attitude  within  two  seconds  after  the  application  of  the 
aileron  displacement.     This  was  caused  by  the  large  roll  angle  that 
developed.     The  resulting  errors  in  the  /&     and     if/      traces   were 
readily  apparent  and  consequently,  these  two  variables  were  disregarded 
in  this  case.     The  only  matching  that  was  attempted  with  aileron  step 
forcing  functions  were  the  amplitude  of  the  initial  peak  of  the  roll  rate 
response,   and  the  neutral  stability  of  the  spiral  mode  as  indicated  by 
the  constant  mean  value  of  roll  rate  after  the  roll  rate  response  had 
reached  its  initial  peak  value.     The  rapid  development  of  the  nose-down 
attitude  when  the  aircraft  was  forced  by  an  aileron  step  function  could 
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have  been  avoided  by  the  use  of  an     initial   condition,    such  as  a  steady 
sideslip.     This  would  have  extended  the  time  during  which  the  airplane 
was  in  the  linear  approximation  range.     Since  the  equations  of  motion 
of  the  analog  are   developed  for   perturbations   from   steady  state,    a 
similar  initial  condition  would  not  be   required  when  forcing  the 
computer . 

After  completion  of  one  series  of  iterations  for  the  values  of  the 
stability  derivatives,  the  procedure  was  repeated  for  all  forcing  func- 
tions,  again  starting   with  the   rudder   pulse.     It  was   noted  that   good 
matching  was   quickly  obtained  for   any  one   of  the  forcing   functions 
taken  individually,  but  a  shift  to  the  other  forcing  functions  offered  a 
relatively  poor  match  in  the  early  stages.     Gradually,  however,   as   a 
result  of  the  experience  obtained,   and  by  reference  to  the  records   of 
the  effects  of  varying  the  individual  derivatives,  new  procedures   for 
establishing  the  value  of  some  of  the  derivatives  became   apparent. 
These  procedures  represented  a  more  definite  method  of  evaluating 
the  corresponding  derivatives  than  did  the  influences  indicated  in 
Tables   ni-a  and  ffi-b,  because  in  each  case  a  particular   response 
characteristic  was  controlled  by  an  individual  derivative.     The  deriva- 
tives and  their  individual  effects  are  as  follows: 

Cn/«         -      The  period  of  the  Dutch  roll  mode  is  controlled  by 
this  derivative.     This  derivative  has  more  than  ten 
times  the  effect  on  the  period  than  does  any  of  the 
other  derivatives. 
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C_  -       C        has  a  controlling  effect  on  the  damping  of  the 

r  r 

Dutch  roll  mode  when  the  analog  is  forced  by  a 
rudder   deflection. 
C  -      Changing  this  derivative  changes  the  effectiveness  of 

by 

the  rudder  in  producing  yaw.     Therefore,     Cn /• 
can  be  established  by  matching   jS    and     (U      ampli- 
tudes  due  to  a  rudder  deflection. 

C\-  -      The  magnitude  of  this  derivative  controls  the  match- 

ing of    0      response  in  two  respects;  phase  and  ampli- 
tude.     Phase  error  in  the  initial    ft      response  of  the 
analog  is  affected  considerably  by  this  derivative. 

C\  .  -       When  the  airplane  is  forced  with  an  aileron  step  func- 

tion,  the   Dutch  roll  mode  is  excited  to  such  a  small 

degree  that  the  initial  peak  of  the   0       response   can 

be  attributed  to  the  rolling  mode  and  the  magnitude 

of    Cw        .      For  a  given    C\        ,     C\  .  can    be 

otu  P  6a, 

determined  by  matching  the  initial  peak  of  the  0 
response.  Determination  of  the  precise  value  of 
Ci  is   explained  below. 

C  q  Just  as    CL      controls  the  damping  of  the    Dutch  roll 

r 

mode,     Cy      controls  the  damping  of  the  rolling  mode. 

With    C\*         adjusted  as  described  above,     Ci       can 
be  determined  by  matching  the    0      response  to   an 
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aileron  doublet,   especially  during  the   latter   half 

of  the   aileron  input. 

Cj   .         -      The   effect  of  this  derivative  is  readily  apparent  in 
Or 

the     0       response  as  an  initial  adverse        0         due 

to  a  rudder  deflection.     Matching  the  initial  form  of 

« 

the     (p       response   will  thus   establish    ^\r 

Cn       )     -       In  the  response  of  the  computer  analog  to  a  rudder 

r.         )  deflection,    Cn_    has  only  a  small  effect.     However, 

cn/    )  P 

)  an  upper  and  lower  limit  can  be  established  for    CL    » 

beyond  which  the  character  of  the       W        trace   is 
noticeably  altered  in  magnitude,    in  the  case  of  an 
aileron  deflection,  the  effects  of    C-      and    Cn  , 
are  so  similar  that  one  can  be  replaced  by  the  other 
to  a  large  extent.     If  an  iteration  is  performed  between 
these  two  derivatives,  however,   maintaining    Cn 
within  the  limits  developed  during  a  rudder  deflection, 
until  the  best  match  is  obtained  in      IP      due  to   an 
aileron  doublet,     Cn      and    Cn*         can  be  determined 
within  reasonable  limits  of  accuracy. 
1  -      Only  a  very  rough  value  of  this  derivative  can  be 

established.     Gross  changes  in    Cj       will  effect  the 
stability  of  the  spiral  mode.      The      0      response  to 
an  aileron  step  function  is  an  indication  of  the   spiral 


r 
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mode    stability:    a  constant   mean  value   of       (p 
indicates   neutral   stability,   a  decreasing  mean  value 
indicates   positive    stability  and  an  increasing  mean 
value  indicates  negative  stability  of  the  spiral  mode. 
Since   the   other  three   derivatives   which  affect  the 
spiral  mode   (  Cn^  ,     Cj,    ,     Cu  )   can  be  determined 
by  the  methods  outlined  above,  the  value  of   Cj      can 
be   adjusted  to  match  the   character   of  the   mean 
response   to   an  aileron  step  function. 
Because   of  the  insensitivity  of  the   responses   to 
changes  in  these   derivatives,  no  refinement  of  the 
theoretically  determined  values  was  possible. 


The   above   relations  are  summarized  in  Table   IV. 


By  the  time  sufficient  information  had  been  accumulated  to 
develop  the  procedures  outlined  above,  the  initial  discrepancies  between 
the  analog  and  the  aircraft  responses,  for  all  four  forcing  functions, 
had  been  appreciably  reduced.     Subsequent  use  of  the  iteration  proced- 
ure,  in  the  case  of  each  forcing  function  and  its  related  responses,   was 
confined  primarily  to  adjustments  in  the  derivatives  which  were  based 
on  the  relationships  of  Table   IV.     In  addition,   when  working  with  the 
responses  to  the  rudder  forcing  functions,   care  was  taken  not  to  make 
large  changes  in  those  derivatives  which  had  been  established  primarily 
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by  matching  responses  to  aileron  forcing   functions.     This  restriction 
was  found  necessary,    since  alteration  of  the  aforementioned  derivatives 
generally  resulted  in  only  minor  improvements  in  the  rudder  cases,  but 
had  an  appreciable  adverse  affect  on  the  matching  of  the  responses   in 
the  aileron  cases.     The  same  restriction  was  applied  when  matching 
responses  to  aileron  forcing  functions  ;    in  this  case,   no  extensive 
changes  were  made  in  those  derivatives  which  had  been  established 
primarily  by  matching  responses  to  rudder  forcing  functions. 

When  a  set  of  derivatives  was  found  such  that  the  responses   of 
the  computer  gave  satisfactory  matching  with  the  flight  data  for   all 
four  of  the  forcing  functions,   and  no  further  improvement  seemed  justi- 
fiable in  relation  to  the  amount  of  time  it  would  involve,  the  final  value 
of  each  derivative  was  read  off  of  the  computer  potentiometers   and 
recorded. 
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RESULTS 

The  final  values  of  the  lateral  stability  derivatives  of  the   NAVION 
as  determined  by  the  method  of  this  report  are  tabulated  and  compared 
with  the  initial  theoretical  values  in  Table   V.     The  matching  of  the 
responses  of  the  final  computer  configuration  to  the  flight  test   data  for 
all  four  types  of  forcing  functions  are  presented  in  Figs.  15-18. 

During  the  iteration  process,   a  better  match  had  been  obtained  to 
the  responses  for  each  of  the  four  forcing  functions  taken  individually, 
but  no  single  set  of  derivatives  gave  as  good  a  match  to  the  responses 
of  all  four  forcing  functions  taken  as  a  group  as  did  the  final  set  of 
derivatives  listed  in  Table   V.     The  major  discrepancies  remaining  in 
the  final  matches  are  the  generally  low  response  of  the  analog  in  the 
/&    variable  for  all  forcing  functions  ;    the  too  large  amplitude  re- 
sponse in  the     tfi       variable  for  rudder  doublet  forcing  functions  ;    and 
the  initial      (1/      response  to  the  aileron  doublet.     It  is  felt  that  had 
more  time  been  available,   further  investigation  into  the  interrelated 
effects  of  the  cross  derivatives    C„      and     Ci       would  have  resulted 
in  nearly  perfect  response  matches  in  the     (If     and     0      variables 
when  using  the  rudder  forcing  functions. 

The  accuracy  of  the  final  value  of  each  derivative  is   directly 
related  to  the  sensitivity  of  the  system  response  to  changes   in  the 

derivative.     The  derivatives    Ci        »     ^1/     '     ^V^   ■     ^n  <      '     ^Irf  ' 

6a,  *r  A  °r  " 
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Ci         and   Gn       are  judged  to  be  very  accurately  determined  by  the 
xp  '  nr 

method.     (See   Appendix    C   concerning   the  accuracy  of  the    specific 

values  for    C\  and   Cn  of  this  report.)    The  values  for  Cl    » 

Cn        and   C  are  susceptible   to  somewhat  larger  error  due    to 

the    smaller    effects   they  produce  on  the  response,    and  also  due  to 
overlapping  effects  in  the    case  of   Cn      and   Cn  •     Finally,    the 

value  of  the  derivatives    Cy         and   Cy  could  not  be  established 

with  satisfactory  accuracy.    Other  types  of  flight  maneuvers  might  be 
advantageous  in  establishing  these  latter  derivatives <> 

The   proposed  method  of  extracting   stability  derivatives   from 
transient  flight   data  was   determined  to  be  completely  feasible.     As 
mentioned  previously,   there  are  inherent  difficulties  in  the  process  of 
extracting  all  of  the   derivatives  from  flight   data  regardless   of  the 
method  used.     For  certain  applications,  however,   no  better  criterion 
could  be  established  for  obtaining  the  values  of  the  stability  deriva- 
tives than  to  match  the  response  of  the  analog  containing  these  deriva- 
tives to  the  response  of  the  actual  aircraft.    In  the  field  of  servomech- 
anisms,  for  example,  this   criterion  would  be   considered  to  be   of 
prime  importance. 

Another  advantage  of  the  method  would  be  its  economy  in  time, 
both  in  obtaining  flight  data  and  in  extracting  the  derivatives  from  that 
data.    In  obtaining  the  flight  data,   any  arbitrary  forcing  function  can  be 
used,   providing  it  offers  sufficient  excitation  of  the  modes  of  aircraft 
motion.      Use   of  the   telemeter    data   collection   system  permitted 
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convenient  monitoring  of  flight   data  during  the  flights,   and  thus  when 
necessary,   instructions  could  be  given  to  the  pilot,  to  improve   the 
responses  obtained.     Normally  this  would  allow  the  recording  of  all 
the  data  for  one  flight  condition  in  a  matter  of  minutes.     The  ease   of 
extracting  the  derivatives  would  tend  to  increase  with  practice.     Un- 
doubtedly,  Table   IV,   which  contains  procedures  for  establishing 
several  of  the  derivatives,   could  be  expanded  with  the  aid  of  more 
time  and  experience  with  the  method.     For  this  reason,   it  is  felt  that 
better  matching  than  is  illustrated  in  this  report  would  definitely  be 
possible.     It  is  important  to  note  that  the  results  of  this  investigation 
were  obtained  for  the  aircraft  in  the  cruise  condition  only<>     Because  of 
the  large  variation  of  certain  stability  derivatives  with  changes   in 
flight  conditions,  further  investigation  of  the  method  throughout  a 
range  of  flight  conditions  is  recommended.     There  is  virtually  no 
chance  of  computational  errors  affecting  the  results,  for  once  the  analoj 
has  been  completed  the  computations  are  performed  by  the  computer. 
Although  the  theoretical  estimates  of  the  derivatives  were  made 
as  accurately  as  possible  in  this  investigation,   it  was  found  that  the 
rate  of  convergence  of  the  iteration  process  was  sufficiently  rapid  in 
the  first  stages  to  suggest  that  very  rough  estimates  of  the  derivative 
values  would  be  sufficient  for  use  as  initial  computer  settings.     The 
time  involved  in  refining  the  rough  estimates  would  more  than  likely 
exceed  the  amount  of  computer  time  necessary  to  accomplish  the  same 
degree  of  refinement. 
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Certain  limitations  of  the  method  are  readily  apparent.      The 
exact  form  of  the  equations  of  motion  must  be  known  in  order  to  con- 
struct the  analog.     The  application  of  the  method  to  more  complex 
aircraft  would  require  further  investigation  and  refinement  ;    for  ex- 
ample,  any  nonlinear  effects  of  the  aircraft  system  would  have  to  be 
duplicated  by  electronic  circuitry  in  the  analog,,     The  accuracy  obtain- 
able for  those  derivatives  to  which  the  response  is  insensitive   is 
undoubtedly  of  low  order,  but  the  possibility  does  exist  of  using  differ 
ent  types  of  flight  tests  to  establish  these  values ,     If  these  tests 
were  compatible  with  the  proposed  data  collection  and  reduction 
method,  the  additional  effort  could  be  minimized,, 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  proposed  method  of  extracting  stability  derivatives   from 
transient  flight   data  was  proved  to  be  feasible  within  the  scope  of  this 
investigation.     Matching  of  responses  between  the  test  aircraft  and  the 
equations  of  motion  containing  the  derivatives  is  an  excellent  criterion 
for  many  applications  of  the  derivatives.     Due  to  the  compatibility  of 
the  data  collection  system  and  the  analog  computer,  the  method  is  rela- 
tively fast  and  free  from  computational  error.     Only  rough  approxima- 
tions to  the  values  of  the  stability  derivatives  are  needed  to  construct 
the  analog.     In  general  the  accuracy  with  which  each  of  the  derivatives 
is  obtained  is  associated  with  the  importance  of  that  stability  derivative 
in  influencing  the  motion  of  the  aircraft  during  the  test.     The  use  of 
four  different  types  of  control  deflections  as  forcing  functions   is   con- 
sidered a  minimum  requirement  for  accurate  determination  of  all 
lateral  derivatives.     Subject  to  these  restrictions,  the  exact  form  of 
each  particular  control  deflection  is  arbitrary.     Application  of  this 
method  should  include  the  preliminary  procedures  described  under 

Development  of  Method  of  Analysis  '    and     Preliminary  Test   of 
Feasibility  of  Method.  " 

Application  of  this  method  to  non-linear  aircraft  would  require 
further  investigation  and  refinement.     Any  non-linear  effects   in  the 
aircraft  system  would  have  to  be  duplicated  by  electronic  circuitry  in 
the  analog. 
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As  a  re  stilt  of  the  investigation,   the  following  recommendations 
are  proposed: 

1.  Investigate  the  method  at  other  flight  conditions  and  with  more 
complex  aircraft. 

2.  Investigate  the  use  of  other  control  deflections  to  obtain 
responses  which  are  more  strongly  influenced  by  the  less 
sensitive  derivatives.     As  an  example,    Cl      might  be  more 
accurately  obtained  from  an  analysis  of  the  response  to  a 
rudder  step  function. 

3.  Apply  an  initial  condition  such  as  a  steady  sideslip,   prior  to 
forcing  the  airplane  with  an  aileron  step  function. 

4.  Augment  the  data  (Table   III)   which  indicates  the  effect(s)   of 
the  individual  derivatives  on  the  response  of  the  system  by 
conducting  a  similar  variation  of  the  derivatives  when  forcing 
the  computer  problem  with  an  aileron  function. 

5.  Improve  the  test  aircraft  aileron  instrumentation  by  recording 
both  left  and  right  aileron  angles  and  converting  to  average 
aileron  angle  by  means  of  analog  computer  scaling  circuits. 
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TABIE  II 


PHXSICAL  CHARACTERISTICS   CF  NAVION 


A.         WIND 

1.     Total  Wing  Area,  Sy  (includes  flaps,  ailerons 
and  19.87  Ft>  covered  by  the  fuselage) 


2o 

Span  by 

3. 

MAC 

1*. 

Angle  of  Incidence 
Root,  ir 

Tip,     it 

5. 

Twist 

Aerodynamic 

Geometric 

6. 

Airfoil  Section 
Root 
Tip 

7o 

Aspect  Ratio,  ARy 

8o 

Taper  Ratio,     A  v 

9. 

Dihedral 

10. 

Root  Chord 

llo 

Tip  Chord 

AHERON  (one  Aileron  oi 

(Beaded  skin,  Frise  Nose  Balance, 
Mo  Trailing  Edge  Extension) 

1.  Area,      Sa 

20  Span,     ba 

3*  Aileron  Deflection,    ^  a 

bo  Boost 

S.  Trim  Tab 


l8bo337  Ft.2 
33o378  Ft. 
68,352  in. 

♦2° 
-10 

2°  31' 
3° 

NACA  hhlSR 
NACA  61*10R 

6.01*1* 

0.$26$ 
7o*° 
7«2    Ft. 
3o92  Ft. 

120°-120O 

Wheel  Throw 
Fixed  bend  tab  on 
right  aileron. 

20161  Ft.2 

61.987  in. 

30°  UP,  20°  DN. 

None 

Fixed  bend  tab  on 
right  aileron  o 


B.         HORIZONTAL  TAIL 

1*     Total  Area 

(Includes  2.368  Ft.2  covered  by  fuselage) 

2,  Span,   by 

3»  Airfoil  Section 

U.  Taper  Ratio,    A  t 

5o  MAC,       /Ct 

6.  Aspect  Ratio,      p^  ^ 

B-lo     HORIZONTAL  STABILIZER 

1.  Area,  Ss 

2.  §h  (with  reference  to  FRP) 

B-2.     ELEVATORS 

(No  Trailing  Edge  Extensions) 


1.  Total  Area,     Se 

2©  Span,     be 

3.  Deflection,     ^  e 

Uc  Trim  Tabs  (32»  Span,  1-1/2 "  Chord) 

5o  Root  Chard 

6.  Tip  Chord 

C.         VHJTICAL  TAIL 

1.     Total  Area,     Sv 

(Includes  2<,£77  Ft.2  blanketed  by 
fuselage  and  excluding  I.I483  Ft.2 
of  dorsal  fin) 

20    Airfoil  Section 
Root 
Tip 


U3.0S1  Ft.2 

13ol72  Fto 
NACA  0012 
.67 

3.31  Ft. 

h.  02 


28.953  Ft.2 
-Uo 

Smooth  Skin 
Flat  sided 
No  trial  bungee 
Balance  spring 

lii.098  Ft.2 
73.582  in. 
30o  UP,  20°  DN. 
30°  UP,  30°  DN. 

1.5  Ft. 

1.0  Ft. 


12.925  Ft.2 


NACA  0013 o 2  Mod. 
NACA  0012 -6I1  » 


C-l.     VERTICAL  STABILIZER 

1.     Area,     Sy  6.873  Ft.2 

'       vs 

(Includes  0.11*27  Ft.2  blanketed  by  fuselage) 
2.06  (With  reference  to  fuselage  0)  2°  Nose  left 


C-2.     RUDDER 


1.  Area,     Sj. 

2.  Rudder  Deflection,   ^  r 

3o  Rudder  Pedal  Throw 

lio  Trim  Tab 


Smooth  Skin 
Rigged  3°  Rt.  to 
Fin  0,    Fixed 

Bend  Tab. 

6.052  Ft.2 
17°  L  -  23°  R 

5.75  in. 
Fixed  Bend  Tab 


D.         PROPELLER 

lo  Activity  Factor,  A.F. 
2.  Diameter,  d 
3«  Pitch  6  i  .75  R 


Hart sell 
100 
86" 

21©  30' 


NOSE  QEAR 

1.  Travel 

2o  Strut  Angle  (To  Vertical) 


♦  30° 
10° 


F.         MISCELLANEOUS 

lo     Length,  overall 
20     Tail  length,     -Avt 


27.25  Ft. 
16.88  Ft. 


0. 


POWER 


1.     Engine 


Continental  E-185  Engine 
185  HP  ©  2300  RPM,  29"  Hg 
M.P.  @  Sea  Level 


TABLE   m-a 

EFFECT   OF  VARIATION  OF  STABILITY   DERIVATIVES 
ON  AIRCRAFT   RESPONSE  W 


PARAMETER 
INCREASE 

PEAK  MAGNITUDE  OF: 

PERIOD 

Tv2 

e 

^ 

f 

S 

decrease 

decrease 

decrease 

- 

strong 
decrease 

Clr 

- 

- 

- 

- 

slight 
decrease 

S 

- 

- 

slight 
decrease 

- 

C1P 

- 

- 

strong 
decrease 

slight 
increase 

slight 
decrease 

% 

- 

- 

strong 
increase 

slight 
decrease 

moderate 
increase 

V 

slight 
decrease 

slight 
decrease 

slight 
decrease 

- 

small 
decrease 

Cn£ 

decrease 

decrease 

decrease 

strong 
decrease 

- 

O/V 

increase 

increase 

increase 

- 

- 

O/V 

- 

- 

- 

- 

- 

°yW 

- 

- 

- 

- 

(1)    Only  rudder  forcing  function  used, 


TABLE  IH-b 


STABILITY   DERIVATIVES   WHICH  INFLUENCE 
SPECIFIC  ASPECTS  OF  AIRCRAFT   RESPONSE 


PEAK  MAGNITUDE  OF: 

PERIOD 

DAMPING 

§ 

v-- 

4> 

s. 

(<V 

r 

<v 

<V 

% 

(%> 

<CV 

(Cnr> 

«v 

<%> 

Cl 

P 

«v 

% 

% 

Derivatives  listed  in  descending  order  of  effect. 

(     )  --    Indicates  a  decrease  in  the  variable  affected. 


TABLE   IV 


PRIMARY   EFFECTS  OF  INDIVIDUAL  STABILITY 
DERIVATIVES  ON  AIRCRAFT  RESPONSE 


CONTROLLING 
DERIVATIVE 


■r 


RESPONSE  CHARACTERISTICS 
Period  of  the  Dutch  roll  mode. 


Cn^ 

Cn  Damping  of  the  Dutch  roll  mode 


Cjj  Magnitude  of      *-p      and     (3      response  to  a  rudder 


"^  deflection. 


Cj  .  Phase  and  amplitude  of  the      <$>      response  to  a 


rudder  deflection. 

Cj,  Amplitude  of  initial      4*      response  to  an  aileron 

^  step  input. 

Ci  Having  established  a  value  for    Ci  .         »     Ci 

can  be  determined  by  the  roll  response   to  an 
aileron  doublet,  particularly  during  the  last  half 
of  the  control  input. 

0j  Adverse  roll  characteristics,   which  are  apparent 

0/u  in  the      <$>      response  immediately  after  a  rudder 

deflection,, 

Cn  An  upper  and  lower   limit  for    Cnp    can  De    . 

^  established  by  an  amplitude  match  of  the     H^ 

responses  to  a  rudder  deflection.      CnD    can  then 
<l,  be  adjusted  within  these  limits  along  with  the 

adjustments  in    Cn  c     to  provide  the  best  match  of 

*f         response  to  an  aileron  doublet. 


(continued) 


TABLE   IV  (continued) 


C^  Response  insensitive   to  this   derivative,   but  a 

qualitive  match  of  the  slope  of  the  mean 
response   to  an  aileron  step  function  is   the   best 
means  available  for  determining    Clr  . 

Cy  Response   too  insensitive   to  these   derivatives 

to  permit  the   determination  of  a   special  effect 


\ 


& 


in  either   case. 


1/ 


TABLE   V 


THEORETICAL  AND   FINAL  EXPERIMENTAL 


LATERAL  STABILITY   DERIVATIVES 


Derivative 

Theoretical 
Value 

Final 

Experimental 

Value 

Cv 

-.592 

-.592 

Cn(3 

.0613 

.101 

Cl0 

-o0814 

-.0595 

% 

-.0510 

-.0705 

% 

-.46 

-.43 

S 

-.0937 

-.0951 

Clr 

.102 

.131 

ClV 

.0152* 

.0069 

%* 

-.0561 

-.0774 

cv 

.111 

.144 

0  Oo 

.133 

.137 

6  Q^ 

-.0087 

-.0031 

*  From  Ref.  9 
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ANALOG    OF    LATERAL  EQUATIONS 
OF   MOTION    OF     THE    NAVION 
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POTENTIOMETER    SETTINGS 
FOR   ANALOG    OF   NAVION 
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APPENDIX   A 

THEORETICAL   AND   EXPERIMENTAL  DETERMINATION 

OF   THE  DYNAMIC  RESPONSE   CHARACTERISTICS 

OF  THE  SIDESLIP   VANE 

General 

The   instrument   system  chosen  for   measurement   of  the    sideslip 
angle  consisted  of  a  light  vane  mounted  on  a  boom  which  extended  forward 
from  the  left  wing  tip.     (See  Fig.  4).     As  has  been  previously  stated,   it 
was  considered  of  prime  importance  that  the  sideslip  measurement  sys- 
tem have  dynamic  response  characteristics  such  that  no  coupling  error 
would  be   introduce  do 

Since  no  data  were  available  for  the  response  characteristics   of 
the  vane  selected,  theoretical  calculations  were  made,  and  experimental 
tests  conducted  in  order  to  determine  if  the  natural  frequency  and  damp- 
ing ratio  of  the  vane   were   within  acceptable   limits. 

Theoretical  Analysis 

(See   Refo  1,  pp  11.  51-11.  55).     Notation  as  per  sketch  below. 

V  -  true  airspeed,  ft/  sec 

V*i  ■  indicated  airspeed,  ft/sec 

M  -  moment  about  pivot 

Cn  »  normal  force   coefficient 

I  •  moment  of  inertia  of  vane 


A-2 


i-pIq 


X 

~  Jf 

Earth  Reference!^ 


Vane  Pivot  A)cis 

Disturbing  Moment 

Aerodynamic 
MorYrio  )   Force,  N) 


Vane  Lifting 
Surface. 

c 

Aerodynamic 

Clr\6rdvoi'5e  ForCe. 


(1)  >-   — -     3  —  (°  However ,   due  to  angular  velocity  of  the 

vane  with  respect  to  the  air  mass,  the  angle  of  attack  at  the  vane, 

A. 


(2)     -x  v    =1  .x  ■+•    y  V  >    +-  r  J 


(3) 


c*v  =    %~(2    +■ 


i* 


&    r 


X. 


(4)        *    =     t  V'^c 


-x. 


AC 


(5)  M  =    t  C-  ^    s  ^v  i^ 

d   _X„ 

(6)  n  -  ^  m  -U+rUo 

combining  (3),    (4),    (5),   &  (6), 


(7)    S  1 


dynamic  pressure 


<&    ^i  +  ■  S  ~a-r*l  -ri 


|\A 


v 


o 


l\ 


I/O 


tvJ 


i  ■> 


A-3 


(8) 


(9) 


where: 


d  *• 


3    =    \     "I 


Ps  T 


lA 


(10)     \<     r 


Jl- 


rvj 


d    * 


Natural  Frequency 
of  Vane 


Damping  Ratio 
of  Vane 


Aerodynamic 
Stiffness  of  Vane 


Conditions  for  Analysis: 

Density  Altitude:      5000  ft 


V 

s 


Jl 


I    : 


120  mph   -  176  ft/sec 
189.6  ft/sec 
.0361  ft2 
.157  ft 

lo72/radian  (See  Ref.   5) 

6.33  x  10        slug  ft      (Determined  experimentally) 


From 


LO 


r~" 


r\ 


From 


176 


U3 


(8)   above, 


.002378     *    .0361      *   .157     x    1.72 


V\ 


\ 


75.5 


U) 


2    » 

Ra 

sec 


6.33     x    10 


12.01    Cycles/ 


sec 


VA 


r 


(9 )   above, 

.002378    v     .8616     *     .  0361-  (,157)3     -    1.72 


8     v    6.33 
.031 


10 


^T 
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Experimental  Analysis 

Since  the  above  theoretical  dynamic  analysis  of  the  response  char- 
acteristics  of  the   sideslip  vane   was  based  on  a  moment  of  inertia  of 
limited  accuracy,   and  neglected  the  effects  of  friction  damping  in  the 
system,  it  was  considered  desirable  to  determine  the  dynamic  response 
of  the  vane  experimentally.     Ref .  1   (p  11.  55)   reported  that  wind  tunnel 
calibration  of  vanes   is   usually  unsatisfactory  because   of  the   poorly- 
damped  vane  response  to  disturbances,   and  the  usual  high  level  of  tur- 
bulence in  wind  tunnels.     For  this  reason,  the  sideslip  vane  system  was 
modified  to  permit  in-flight  dynamic  testing  of  the  vane  when  installed 
on  the  airplane.     A  spring  and  plunger  mechanism,   actuated  by  a  sole- 
noid,  was  installed  on  the  boom,   as  indicated  in  Fig.  4.     When  in  flight, 
operating  at  the   desired  test  configuration,  the  vane  was   deflected  by 
sideslipping  the   aircraft,   and  locked  from  the   cockpit  by  electrically 
extending  the  plunger.     (See     Solenoid  Switch,      Fig.   3.)    Thus,   when  the 
solenoid  was  reenergized,  the  vane  was  released,   and  responded  as  if 
to  a  step  function.     Fig.  A-l  shows  the  transient  response  of  the  vane  to 
such  a  step  input.     Calculations   for  the  experimental  determination  of 
the  natural  frequency  and  damping  ratio  of  the  vane  from  the  transient 

response  of  Fig.  A-l  are  as  follows: 
Conditions  for  Test 

Density  Altitude:       5000  ft 
Vi  :        120  mph    ~   176  ft/sec 
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In  Ref.   6,   it  is  shown  that  the  undamped  natural  frequency  and  the  damp- 
ing ratio  may  be  determined  from  transient  data  as  follows: 


^_ 


where 


Sv  J?         X2 


■     logarithmic   decrement 


From  Fig.  A-l, 

"X"  -     .065   sec,    Xx    -=    4.0,      X2     »   3.1 


U 


*\ 


(.2546)2  +         IT2 
\T  .065/2 

r         i£n_ 


b*^ 


.2546 
3.15 


97.0  Ra/seCo 

15.44   Cycles/ sec 
.081 


Comparison  of  Results 
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fn,  Natural 
Frequency 
cps 


•^  ,   Damping 
Ratio 


Theoretical 
Analysis 


12.01 


0.031 


Flight 
Test 


15.44 


0.081 


Concluding  Remarks 

(1)  The  natural  frequency  of  the  sideslip  vane  has  been  shown  to 
be  of  the  order  of  16  cycles  per  second,   and  the  damping  ratio 
approximately  0.08. 

(2)  The  ratio  of  the  natural  frequencies  of  the  sideslip  vane,   and 
the  airplane  in  the   lateral  oscillation  is  approximately  40. 

(3)  In  view  of  (1)  and  (2)   above,  no  dynamic  correction  to  sideslip 
vane  data  was  required. 


o 

Ld 


LU 
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APPENDIX  B 

DETERMINATION    OF    LATERAL    STABILITY    DERIVATIVES 
FOR  THE  NORTH  AMERICAN  NAVION  FROM  THEORETICAL 
FORMULAE,   AND   FROM   WIND  TUNNEL  AND  FLIGHT  TESTS 
FROM  GEOMETRICALLY  SIMILAR  AIRPLANES 


The   determination  of  stability  derivatives  for   an  aircraft  from 
theoretical  and  experimental  data  is  a  problem  which  has  received  wide- 
spread attention.,     Numerous  publications,   outlining  a  considerable  num- 
ber of  methods  for  calculating  or  estimating  derivatives,   and  presenting 
the   results   of  extensive  test  programs,  have   appeared  over  the   past 
twenty  years.     (See  bibliography  to  Ref.  7).     However,  in  spite  of  such 
a  wealth  of  information  on  methods  and  test   data,   estimation  methods 
are  often  found  in  practice  to  yield  only  fairly  accurate  values,   suitable 
for  making  only  first  approximations   to  dynamic  stability  characteris- 
tics.    In  addition,  the  application  of  several  theoretical  methods  to  the 
same  airplane  often  results  in  wide   discrepancies  between  some   of  the 
calculated  derivatives. 

In  general,  the  most  reliable  methods  for  calculation  of  lateral 
stability  derivatives  depend  on  wind  tunnel  data  for  accurate  values  of 
the  sideslip  derivatives,     v-^  ^   ^       nfi      ^         J?  A  •     This   is  par- 

ticularly true  for  propeller  driven  aircraft,  where  the  propeller  lateral 
force  and  slipstream  effects  are  important  contributions  to  the  sideslip 
derivatives.     No   satisfactory  purely  theoretical  method  has   yet  been 
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developed  for  accurately  estimating  these  derivatives  for  a  complete  air- 
plane primarily  because  of  large  interference  effects  between  the  various 
airplane   components . 

The  remarks   above  concerning  wide   discrepancies  were   found  to 
be  particularly  pertinent  to  the  calculation  of  the  lateral  stability  deriva- 
tives of  the  North  American  NAVION.     The  derivatives  for  this  airplane, 
including  the   sideslip  derivatives,  were   calculated  by  five   well  known 
methods,   and  the  results  compared  with  test  data,  both  for  the  NAVION, 
and  also  for  several  geometrically  similar  airplanes.     The  results,  pre- 
sented in  Table  B-:ILJ  showed  a  variation  of  at  least  100%  between  the  high 
and  low  values  in  the  case  of  all  derivatives  with  the  exception  of   ^  n  » 

which  varied  by  approximately  30%. 

It  is  to  be  noted  in  Table  B-II  that  the  method  of  Ref.   5   gave  the 
best  correspondence   between  theoretically  calculated  derivatives,    and 
the    "average"  derivatives   of  the   several  methods   and  tests   results. 
Therefore,  the  detailed  calculations   for  the  lateral   stability  derivatives 
of  the   NAVION  by  the  method  of  Ref.   5  are  presented  in  the  following 
pages. 
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GENERAL   PARAMETERS   FOR   NORTH  AMERICAN 

NAVION  N91566 

In  addition  to  the  physical   characteristics   of  the   NAVION  pre- 
sented in  Table   II  to  the   basic   report,  the   following  parameters, 
determined  from  the   assumed  test   conditions,   are  required  for    cal- 
culation of  stability  derivatives. 


W    . 

V    . 
CL 


.  Gross  weight  during  test,   lb» o  .  .  .  2770 

.  Density  altitude,   ft  ............ 5000 

.  Air  density,   slug/  ft3  .  . . .002049 

,  Indicated  airspeed,   mph 1*0 

.  True  airspeed,   ft/sec ,  189.  6 

c  Total  lift  coefficient 0 .  407 

.  Tail  efficiency  factor  (estimated)  .......  0.90 


Unless  otherwise  specified,   all  further  references  in  this  appendix 
are  to  Ref.   5  to  the  basic   report. 
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TABLE  B-II  (a) 

THEORETICAL  DETERMINATION  OF  THE  LATERAL  STABILITY 
DERIVATIVES  CF  NAVION  N91566 


1    - 

1 

1  - 

1 

-— 



i  Reference 

J 

i 

(5) 

Per kin* 

&  Hage 

-.592* 

(8) 
USNPGS 

Notes 

-.592 

1 

TR 
635 

(12) 
TR 
589 

(13) 

TN 

1581 

AVO. 
VALDE 

C^ 

• 

1 

-.632 

1 

-.612 

1 

I         C 

|      .0613 

1 

I 

.0708 

1 

- 

.01409 

- 

.057 

C*p      j    -.0816 

1 

-.1236 

i 

-.0981 

-oll37 

-o0991 

j 

-ol07 

1 
-.195 

-ol*6 

-.16 

! 

!      -.301* 

-.l»ii 

Cf          j      .102 

.110 

- 

.102 

.0978 

.103 

On-       !    -.Q51 
*       ! 

-002I43 

- 

-.01*31 

-.0221* 

-o035 

Cn^         -.0937      : 

-.091*5 

- 

! 
-.0953 

-.095 

C            ' 

.111 

- 

- 

j 

«. 

.111 

! 
^    Si 

i 

i 
-.0087 

- 

* 

\ 

" 

! 

-.0087 

c»v 

-.0561 

- 

- 

1 

_ 

-.0561 

c'.,  i 

.133 

" 

"* 

- 

- 

.133 

*  From  Ref , 

.  8. 

TABIE  B  II  (b) 

COMPARISON  OF  LATERAL  STABILITY  DERIVATIVES 

FOR  AIRCRAFT  GEOMETRICALLY  SIMILAR 

TO  NAVION  N91566 
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REFERENCE 

(9) 

P.R.* 

232 

(13) 

n 

2195 

(HO 

TN 
1327 

(15) 

RM 

L53K09 

(16) 

RM 

A52J06 

AVG 

VALUE 

cL 

.390 

o33li 

.Uo 

.33 

.1*12 

sVs„ 

.068 
.506 

.088 
.1*60 

.185 
.1*1 

.131 
.1*56 

.10 
.1*51 

6.0b 

6.2$ 

5.9 

5.89 

6.25 

A 

.527 

.503 

.500 

.$2 

.576 

.  ^M- 

7.5 

6.0 

60O 

3.0 

2.0 

A-Jwj- 

0 

0 

0 

0 

0 

STABILITX  DKRIVA 

TIVES 

r 

-.595 

-.802 

-.573 

-.515 

-.638 

-.621* 

^ — n  n> 

.11*1 

.132 

.0975 

0 111*6 

.109 

.120 

c^ 

-.0822 

-.08 

-.0688 

-.0573 

-.078 

-.073 

C*f 

-.1*3 

-.39 

-.1*76 

-.1*32 

C-Aj^ 

.03 

.06 

.U4O 

.073 

v_  n  ^q 

-.007 

-.02 

-o020 

-.018 

^/v 

-•160 

-.13 

-.138 

-.11*3 

Cvw 

.203 

.203 

A  in., 

.0166 

0OI66 

-.107 

-.107 

C*W 

.101 

.101 

c 

5  GL- 

-.0035 

-.0035 

♦  For  Navion  5113 
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NOTE:    Ref.   5  does  not  present  a  theoretical  method  for  calculation  of 
MHi  .      Therefore,      ^^    was  calculated  from  the  method  of 
Ref.   8. 

^  c:  where 

L    •      Fuselage   length,   Ft, 27.25 

C.         „      -.12  (33.38)  (27.25)  _ 

Y  184-2  -^=- 


c 

Total  directional  stability,  ^a    ,  is  composed  of  contributions 
from  the  wing,  fuselage,  propeller,   vertical  tail,   and  interference  cor- 
rections  due  to  wing  configuration  and  to  sidewash  from  the  fuselage. 
In  equation  form, 

airplane 


'^^^wCSW^vt-  *-cv  *»CV 
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(a) 


c, 


e 


w 


-.0006  \A\        x   57.3  =  0 


(b)    A,,Ck\/2>      ■       0  for  low  wing  aircraft    (see  p.   320) 


e 


(c)  c^ 


/2>  fuselage:  An  empirical  relationship  for  this  contribu- 
tion, which  is  particularly  applicable  to  single  engined  aircraft, 
is  as  follows: 


W  ••'%(*)(*)©'© 


h 


where 


Ss    «   projected  side  area  of  fuselage  in  Ft 
Lf    «  fuselage   length 

Lj,   L2,    W|,    W2    from  accompanying  sketch 
Kp,      from  Fig.  8-4,  p.   320 
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For  the  Navion,   the  above  quantities  are  as  follows: 

SB         ft2 96.22 

Lf'       ft •    •   .   .   .   .   .   .    .        27.25 

Ll*      ft • 4.25 

L2'      ft    ...... 2.61 

Df'       ft 5.51 

dl»       ft • .         4.99 

^max»  ft 4.40 

Wl"      ft • 4.  14 

W2'     ft 1.76 

&2»       nosetoc.g.,  ft 8  05 

K^f    (from  Fig.   8-4) „uo 

therefore, 

U,  -      -(.96)  (.U0)/96. 22^/27.25^/4.25X7,   76 

f  fuse  >{  184#  2  )  (  33758 /lT?rj  (*TT 


^/3  -    -.0434 

(      fuse 


^C    proP        from  pp.   321, 


0  -IT  T}  X,       — ^  M 

4  SM  k 


w     '  — ; 


D      =     Propeller  diameter,  ft 7oq 

N      »     Number  of  propellers „        i 

Jf^,    -     Distance  from  e.g.  to  propeller  plane,  ft;    7.90 
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J  ^y  «\ 


H/^  V      L  *"  ^Tc^-o  where,  from  p.   321. 


k     ■  1.  5    for  full  power, 
s 

is 


c\  CN  > 


-     .  00165  for  two  bladed  propeller 


For  Navion  N91566  at  full  power, 
SVP  -    185.       ^f  -    .632      (Ref.  9) 

"Y  -    550^    SHT  550  x  .632  x  185 

o   ,  U'        "  .002049  x  (189.  6)3  x  (7)2 

I         ■                          0940 
'  C     full  power       - — - 

For  cruise  power, 

S>K       .  130,  '^1^  .    .78 

-r-  550  x  .78  x  130  0gl5 

X  C      '  .  002049  x  (189.6)3  x  (.  7)2 

k  cruise  power     -     j£  cruise  Power     K   (1.5) 

I  _     full  power 


.0815  „ 

x  1.  5  1.  3 


cruise  power         .0940 


51, 

cruise  power 


JC  \  "  l*3  (-00165)  ■  .  °°2i4 

d  «f 

P  -Tl"(7)2  (7.90)  (.00214)  (1)  (57.3) 

^-n/>>      prop"  4  x  184.2  x  33038 

p  .      -.00610 

^-AA      prop 
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<e>      ^nP  vertical  tail 

Q.^,  the  slope  of  the  side  force  curve  for  the  vertical  tail,  is  a 
function  of  the  effective  aspect  ratio  of  the  vertical  tail,  /^\  . 
From  p.   325, 

fO         ■    1.55       v"t        where    Svt    is  the  nominal  vertical  tail  area 

as  indicated  on  p.   325. 

Q  .   I- 55  x  (4. 05)2      _    2   ^ 

rV\  eNavion  10.89 

From  Fig.   8-8,  p.  324, 
(Xvt    -    -048/deg    .    2.75/radian 
J/  ^    ■     distance  from  airplane  e.g.  to  aerodynamic  center  of 

the  vertical  tail. 

As  measured  from  a  scale  drawing  for  appropriate  e.g., 

J?    .  -   16.88  Ft. 

^  vt^avion 

vt    ■       —p-     '       assumed      ■       .90 
To 

When  sidewash  effects  are  included  as  a  part  of  the  tail  efficiency 
factor  instead  of  as  an  increment,        /-^>  x  C  y-\  ^     ,  the  total 
efficiency  factor  can  exceed  unity,   and  is  of  the  form: 

A  CT  \   rr\  See  for  example.   Fig.  21a  of  Ref.  10. 


n 


Cft^  2.75  (12.925)  (16.88) 

n  p    *  (184.2)  (33.38)    l     ' 
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CO    l\  z  (~r\(^.  »     sidewash  or  interference  flow  correction  due  to 
wing -fuselage  combination. 
From  Fig.   8-7,  p.   325,   with  Df  and  d,  from  part  (  c )   above: 

<H/Df      -    J±2i      «  .904 

b-51 

l\zLVl(1      -    .0004x57.3     .    .0229 
therefore, 

r  \  wing  \  fus  V  \  prop 

-     0      -.0434     <0      -.0061    +.0877    +.0229 
A  /g     .      .  0613 


C-e<o 


^p>  airplane  \   wing  I  vert,  tail  V 

where 


Cx    \^'j)  -  is  a  correction  for  wing-fuselage  interference  due  to  wing 
\ 

height. 

/\  ^  vj  C  A   is  a  correction  for  wing-vertical  tail  interference  as  a 

function  of  wing  height. 

n 

The  effects  of  power  on    \  p   _    are  negligible  for  the   "clean" 

T 

condition. 
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(a)  Q  r\ 

s  P  wing:  (See   Fig.  9-3,  p.   344) 

For      rv\      .      6  Extent  of  dihedral 

A      -    .5  -  1.0, 

L_  wing,   due  to  P      "     *'      /radian 

r 

&  C^:      tipshape      "      °»     from  Fig°  9-5,  p.   345. 
^  /\  ,  C  y  :     From  E*!*  9-6,  p.  347. 

=      c  0008x57. 3    „      .0458 
s  '*  -V/Q  low  wing 


(c)  Co 

"K  p  vertical  tail 
C?^  vertical  taU.-a^^i^ 

where    Z^    »    distance  from    X    axis  to  aerodynamic  center 

of  vertical  tail. 

Z^.  for  Navion      «     4.10  Ft. 

Co  -    -2.75    (12-925)    -LidSL     (090) 

^A    vt  C184-2)      P3.38)      <°-VU' 

^Pvt     '       -0214/ra 
(d)  A2  C^^   :      From  Eqn.  9-11,  p.   346: 

^  2Q  "      -.00016  x  57.3 


low  wing 


(\   zQr,  -     -■  00916 

"*/3   l°w  wing 
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Airplane  *  p  wing         '^  vt  *(*  (* 

-      -.0968      -.0214     +  .0458      -.00916 


J>                         •     --08l6/ra 
Vo    airplane /  ra 


■jp    >     Damping  in  Roll 


See   Fig.  9-14,  p.   351. 

For     X    -  0.5.       A^    -    6.04, 

c 

^vo  .     -.46 

0      Navion         ■ 

Referring  to  pp.  428-429 

C0 


A /i  ^  Roll  Due  to   Yaw 


By  strip  integration,   and  assuming  elliptical  lift  distribution, 
Co  CL  .407 


*T  j  Yaw  Due  to  Roll 


^  ^    "    "~g-»  also  by  strip  integration  and  for 


elliptical  lift  distribution. 

^Pvo  ..407 

*  *        — o~    "       -»051 

Navion  °  ■ 
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c 

'  n  so  >    Damping   in  Yaw 


CV     .vCdw-3/J..^1    f^l) 


«„,   -     _^=-  "^vt-g-     ^__;      n^ 


Dw  :  From  p.  95, 

2 


CD  .       CD  ^L_ 


CDW  "        •008 

ffo 

(^       -       6.04 

^       -         .78 

r    z 
CDW  -         .008       +     .07  °L 

For    Cl   ■       .407,      CDw      »      .0196 

Cn  .       -.0196      -2M2.7.)y(^.925)/l6.88\2 

Vb  4  (184.2)  \33. 38/ 

-.0049     -  .0888 


Cv- 


nA, 
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^n  ^ 


where        j       ■     rudder  effectiveness, 
From  Fig.   5-33,  p.  250 

srudder      „        6.052       _        ,,q 


Svt 

12.925 

T  - 

.64 

^  n  s. 

0  ru> 

»    -2.75    (.64 

c 

(184.2)     V33.38 
-.0561 


^  5>~ 


^      Si  "  ^-V-N     c 

^        A  vt        * 

r^Navion       From  Ref.  9      .       .155 
C^  »    -.155    (-.0561)     »    .0087 

C 


V 


t. 


(0.9) 


CV  --£rc"w     "1§^<-^> 


^ 


111 
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^  *&  <x.      ■  "    ^r    X    ^ 


J^  ■    aileron  effectiveness,   a  function  of  aileron -wing 
chord  ratios. 

C    \ 

^M      ifl  a  function  of  the  location  of  the  inboard  and 

"T<v  J 

outboard  ends  of  the  aileron  in  percent  of 

semispan.      (See   Fig.  9-3,  p.   344.) 
Aileron  inboard  end  7.74 

JL^±-        «,   .42 
Aileron  outboard  end  11.4 

Wi -  TT9l    "    -955 


^     .   .70 


Net     S^*2-  -     -™  "   -42 

From  Fig.  9-15,  p.   358, 


.28 


For    : — £-       .    .284,     and  for 
Frise  ailerons,  ^  «       .475 


a. 


c 

^^c  -     .475    x    .28 

V.    B  -       .133 
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Jx  '      ^H 


Jx  2       X        *        kx     -    -J*- 

r  b  -      m 


JB  -       2       kB     ''    .       k2         .j| 

8        -  b  ■     =     m 


From  Ref. 

3, 

for  the   Navion, 

t 

s 

1060   slug -ft2 

I» 

a 

3552   slug -ft2 

W 

= 

2770  lb. 

Jx      -      21 

X 

2  x  1060 

7      -      .0222 
^2  2770/32 . 2      (33.  38)2  ===== 


mb 


mb 


2 


2  x  3552 
2770/32>2    (33. 38)2 


C-l 


APPENDIX  C 
CORRECTIONS  TO  RECORDED  AILERON  DATA 

During  preliminary  matching  of  responses  of  the  test  aircraft 
and  the  computer  to  aileron  forcing  functions,  it  was  determined  that 
the  recorded  values  of  aileron  position  angle  would  require  correction 
for  errors  introduced  as  a  result  of  a  combination  of  aileron  cable 
stretch  and  a  non-linear  hinge  moment  curve.     In  order  to  quantita- 
tively establish  the  extent  of  the  required  correction,  the   computer 
problem  was  forced  by  an  aileron  step  input  to  the  left  and  the   roll 
rate  response  matched  to  the  flight  data  by  adjusting    ^1  r        ; 

On, 

a  similar  aileron  input  to  the  right  was  then  introduced,   and  was 
observed  to  cause  a  computer  response  in  roll  rate  that  was  smaller 
than  the  flight  data  by  a  factor  of  0.78.     This  factor  of  0.78  was 
proved  to  remain  constant  over  a  wide  range  of  values   of     Cl_  , 
the  one  stability  derivative  that  would  have  a  marked  effect  on  the 
roll  rate  due  to  aileron  deflection.     Similarly,  the  roll  response  to 
an  aileron  doublet  could  be  brought  into  perfect  match  with  the  flight 
data  during  the  initial  roll  to  the  left,  but  no  reasonable  value   of 
Ci       would  produce  as  high  a  roll  rate  to  the  right  as  was  recorded 
in  the  flight  data.     (Values  of    Cj      as  low  as    0.28  were  investigated.) 

At  this  point,   equation  (9-33)   of  Ref.   5  was  rewritten  in  ratio 
form,  and  applied  to  responses  of  actual  flight  data  of  aileron  step 
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functions  to  the  right  and  left,   as  follows: 


0  left        _      j    Sep  left 
(p   right  ^      hcL  right 


The    k    factor   was  inserted  in  the   equation  knowing  that  its   value 
should  be   unity.     However,   upon  solving  for    k    it  was   again  deter- 
mined to  be  the  factor   0.78. 

This  analysis  offered  sufficient  proof  that  an  instrumentation 
error  had  been  included  in  the  data.     The   gyro  measuring  the    roll 
rate   response  was  known  to  be  highly  accurate,  but  aileron  deflec- 
tion was  measured  directly  at  the  right  aileron  and  the  calibration 
curve  was  obtained  on  the  ground  under  a  no-load  condition.     Aileron 
cable  stretch  was  known  to  exist  to  some  degree  in  the  aircraft  used, 
and  for  this  reason,   aileron  deflection  was  measured  directly  at  the 
right  aileron.     Instead  of  alleviating  the  instrumentation  error,  how- 
ever, this  merely  caused  it  to  be  additive  during  left  roll  and  sub- 
tractive  during  right  roll.     This  condition  can  be  demonstrated  by 
first  considering  the  non-linear  hinge  moment  characteristic  of  a 
modified  frise  aileron: 


+  a 
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It  should  be  noted  that  the  change  in  hinge  moment  is  nearly  pro- 
portional to  the  deflection  when  the  aileron  is  displaced  down.     There- 
fore,  due  to  cable  stretch,   the  down-going  aileron  would  be  attenuated 
by  a  constant  factor  in  comparison  to  the  no-load  calibration.     The  up- 
going  aileron  has  a  relatively  constant  hinge  moment  over  the   range 
used  in  this  investigation  and  would  therefore  parallel  the  no-load  cali- 
bration curve.     The  resulting  error  is  illustrated  by  a  plot  of  control 
deflection  vs.   aileron  deflection: 


(L£FT) 


(em  nr) 


no-  load 

CALIB#AT/0AL 

ACTUAL 

UNDFX    LOAD 


When  the  ailerons  are  deflected  to  the  point  where  the  right  aileron  is 
at  position  A,  the  signal  output  is  proportional  to  the  length  1-3,   but 
the  actual  deflection  is  the  shorter  length  1  -Z.     Correspondingly,    if 
the  controls  are  reversed  to  the  point  where  the  right  aileron  is   at 
point  B,   the  signal  output  is  proportional  to  the  length  6-7,    but  the 
true  position  is  then  the  larger  value   4-5. 
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Qualitatively,   this  discussion  would  hold  true  whenever  the  slope 
of  the  hinge  moment  curve  differs  between  the  up-going  aileron  and  the 
down-going  aileron.     Although  no  hinge-moment  data  was  available  for 
the   Navion,   the   ailerons   are   of  the   modified  frise  type,   which  are 
known  to  have  this  characteristic. 

Having  qualitatively  explained  the   error  in  terms   of  non-linear 
hinge   moments,   and  having  quantitatively  found  the   magnitude   of  the 
error  both  from  flight  data  analysis  and  from  the  computer  responses 
to  aileron  step  functions,   it  then  became  necessary  to  remove  the 
error  from  the  signal  in  order  to  continue  the  computer  analysis. 
Although  the  ailerons  had  a  differential  gearing  characteristic  for  large 
deflections  which  would  affect  the  apportionment  of  the  error  between 
the  right  and  left  signals,  the  differential  effect  was  less  than  ten  per- 
cent in  the  small  range  of  deflections  used,   and  thus  the  apportionment 
of  the  error  could  be  made  on  a  50-50  basis  with  reasonable  assurance. 
(This  concept  could  be  verified  by  making  reference  to  the  plot  of  con- 
trol deflection  vs.  aileron  deflection.     With  no  differential  in  the 
aileron  system,  the  percent  of  error  for  a  given  control  deflection  to 
the  left  is  equal  to  the  percent  of  error  for  a  corresponding  control 
deflection  to  the  right,   but  the  signs  of  the  errors  are  opposite.     Adding 
a  differential  gearing  effect  changes  this  equality  between  the   error 
signals,   but  not  to  a  marked  degree. ) 

To  effect  the  change  in  the  aileron  signal,   circuitry  was  provided 
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to  divide  the  signal  into  two  paths,   one  path  for  the  signal  during  right 
roll,   and  another  path  for  the  signal  during  left  roll.     The  right   roll 
signal  was  then  attenuated  by  the  previously  established  factor  of  0.78„ 
This  circuitry  is  shown  in   Fig.   9. 

With  this  change  the  character  of  the  signal  had  been  corrected, 
but  the  magnitude  of  the  entire  signal  was  known  to  be  low  because  the 
correction  had  all  been  placed  in  one  direction  as  an  attenuation  rather 
than  being  apportioned  evenly  between  the  two  directions  of  roll  as   an 
attenuation  of  one  and  as  amplification  of  the  other .     The  amount  by 
which  the  signal  was  low  was  the  ratio  of  0.  7 8/0 „  89,   or   (0.877),   for 
dividing  the  signal  by  that  factor  would  halve  the  previous  attenuation 
in  the  right  signal  and  amplify  the  left  signal  correspondingly.      Thus 
the  correction  would  be  apportioned  equally  to  the  two  signals: 


&%         16% 
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Instead  of  applying  this  second  correction  factor  to  the  signal  which 
would  have  required  additional  changes  in  the  circuitry,  the  investiga- 
tion was  carried  to  its  conclusion  and  the  correction  factor  was  then 
inversely  applied  to  the  value  of  the  aileron  derivatives,  the  net  effect 
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of  the  two   alternative   methods   of  correction  being   identical   in  the 
final  result. 

In  conclusion,   it  should  be  stated  that  the  corrections  made  to 
the  signal  were  justified  by  the  viewpoints  brought  out  in  the   above 
discussion.     However,  to  state  that  the  error  was  completely  elimin- 
ated would  seem  overly  optimistic,   especially  with  respect  to  the 
apportionment  between  attenuation  of  the  signal  in  one  direction  and 
amplification  of  the  signal  in  the  other  direction.     It  is  strongly  felt, 
however,  that  any  error  remaining  in  the  aileron  signal  had  no  effect 
on  the  value  of  any  of  the  derivatives  other  than  the  aileron  derivatives, 
Ci  and    C_  •      The  magnitude  of  the  error  produced  in  these 

two  derivatives  should  not  exceed  five  to  eight  percent. 

A  more  suitable  method  of  instrumentation  would  have  been  to 
include  a  transducer  at  each  aileron,   and  sum  the  two  aileron  angle 
signals  and  compute  the  average  angle  as  a  function  of  time  on  the 
analog  computer. 
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